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Attenuation  of  infrared  radiation  in  the  troposphere  is 
dominated  by  water  vapor  absorption.  Most  past  work  on  the  spectro¬ 
scopy  of  water  vapor  has  been  on  the  analysis  of  the  rotational 
and  vibrational  bands.  As  a  result  a  thorough  listing  of  line 
positions,  line  strengths  and  halfwidths  exists  today.  Because 
of  the  weak  "continuum"  absorption  in  the  water  vapor  windows 
centered  at  10  and  4  microns,  efforts  to  measure  and  model  the 
continuum  have  not  been  as  successful  as  efforts  in  tho  analytic 
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of  the  bands.  Yet  a  precise  understanding  of  continuum  absorption 
is  important  for  long  path  energy  transmission.  For  this  reason 
the  study  of  the  water  vapor  windows  has  had  a  long  and  speculative 

history,,  ...  - - - 

^4/ 

The  purpose  of  this  study  is  to  demonstrate  the  importance 
of  far  wing  phenomena  in  characterizing  H~0  continuum  absorption. 

A  total  line  shape  for  water  vapor-nitrogen  interactions  valid 
under  tropospheric  conditions  is  derived.  The  model  is  tested 
using  a  set  of  experimental  room  temperature  H?0  continuum  measure¬ 
ments  of  high  quality.  Using  this  data  base  parameters  of  the 
far  wing  component  of  the  total  line  shape  are  determ ined,f£om 
near  band  experimental  data.  Grating  spectrometer  measurements 
from  300  to  650  cm"  are  used  to  determine  unknown  far  wing  param¬ 
eters  of  the  pure  rotational  band  of  H~0.  CO  and  HF  laser  measure¬ 
ments  taken  in  the  5  micron  and  3  micron  regions  are  used  to 
determine  the  far  wing  parameters  of  the  v?  and  fundamental 

bands,  respectively. 

The  total  line  shape  model  is  applied  to  the  millimeter, 

10  pm,  and  4  pm  absorption  windows  with  encouraging  success. 

A  significant  increase  in  the  self-broadening  ability  of  H?0  over 
N-  is  predicted  in  the  far  wing.  This  allows  the  proper  modeling 
of  the  absorption  coefficient  versus  H^O  partial  pressure  dependence 
in  all  window  regions.  A  negative  temperature  dependence  is  predic¬ 
ted  by  the  model  in  the  continuum.  The  observed  rate  of  the  tempera¬ 
ture  decrease  is  not  predicted  by  the  model,  however  the  failings 
are  related  to  approximations  made  on  the  interaction  potentials 
and  perturbation  expansion  of  the  Hamiltonians.  Although  the 
total  line  shape  model  has  limitations,  it  does  demonstrate  the 
importance  of  considering  far  wings  of  absorption  lines  in  continuum 
absorption. 
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WATER 


CHAPTER  I 
WATER  EVERYWHERE 


Rhyme  of  the  Ancient  Mariner 
-  Samuel  Taylor  Coleridge 


The  importance  of  water  is  suggested  by  its  considerable 
abundance  on  the  surface  of  the  earth.  Its  physical  properties 
have  been  studied  for  a  long  time  and  it  has  proven  to  be  a  molecule 
of  great  complexity  which  gives  water  special  and  unique  properties. 

One  of  the  largely  unsolved  problems  concerning  atmospheric 
water  vapor  is  its  effect  on  infrared  (ir)  propagation  through 
the  atmospheric  windows.  The  propagation  of  radiation  is  affected 
by  three  major  phenomena: absorption,  scattering  and  turbulence. 

In  the  atmosphere  water  is  primarily  in  its  vapor  phase;  however 
water  vapor  clusters  can  form  on  particulates  (i.e.,  clouds,  fog) 
or  ionized  molecules.  Attenuation  of  ir  radiation  by  molecular 
water  vapor  is  largely  due  to  absorption.  Water  vapor  clusters 
contribute  to  both  absorption  and  scattering.  Local  heating  of 
the  atmosphere  caused  by  water  vapor  absorption  can  contribute 
to  turbulence  or  beam  dispersion  effects.  An  ideal  approach  to 
solving  the  propagation  problem  is  to  separate  each  phenomenon 
for  independent  study.  To  some  extent  this  can  be  achieved.  Under 
larboratory  control,  turbulence  and  particulate  scattering  can 
be  minimized.  However  pure  water  vapor  samples  always  contain 
some  dimer  and  cluster  formations  [1]. 

Infrared  absorption  by  the  water  vapor  molecule  is  the  dominant 
mechanism  of  atmospheric  ir  attenuation.  Figure  1  shows  the  low 
resolution  infrared  transmittance  of  the  atmosphere  and  demonstrates 
the  importance  of  water  vapor  over  other  atmospheric  constituents. 
Water  vapor  cluster  attenuation  caused  by  absorption  and  scattering 
is  much  smaller  than  water  vapor  molecular  absorption  in  the  strong 
absorption  band  regions. 

This  study  will  concentrate  on  the  water  vapor  window  absorp¬ 
tion  problem  under  tropospheric  conditions.  The  main  rotational 
and  vibrational  bands  have  been  extensively  characterized  by  Benedict, 
Calfee,  Camy-Peyret  and  Flaud  [2,3,41.  This  work  has  resulted 
in  a  compendium  of  absorption  line  parameters  maintained  by  AFGL  [5]. 
This  compendium  represents  a  significant  contribution  to  absorption 
calculations  and  is  used  in  this  study. 
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lor  general  electro-optical  system  applications  the  strongly 
absorbing  regions  are  not  of  interest;  instead  the  regions  of 
small  attenuation  are  important.  The  infrared  windows  of  the 
troposphere  occur  between  13  to  8  pm,  5  to  4.5  pm  and  4.1  to  3.2  pm 
bounded  by  strong  H~0  or  CO.,  bands.  The  weak  absorption  which 
does  occur  in  the  window  regions  can  be  described  as  arising  from 
two  sources;  local  line  and  continuum  absorption. 

Weak  absorption  bands  of  C0o,  HDO  along  with  other  H~0  absorp¬ 
tion  lines  in  the  window  regions  comprise  the  local  line  contribution. 
The  continuum  contributes  a  slowly  varying  frequency  dependence 
to  the  absorption.  Elsasser  (7),  in  1938,  recognized  the  existence 
of  a  continuum  in  the  13  to  8  pm  window  region  which  he  attributed 
to  the  far  wings  of  the  strong  nearby  rotational  and  \>-  vibrational 
bands  of  H„0.  Further  verification  of  this  nonlocal  line  absorption 
feature  was  provided  by  Yates  and  Taylor  [8]  who  studied  infrared 
attenuation  along  horizontal  paths  at  sea  level.  Solar  spectra 
studies  also  indicated  continuum  absorption  in  the  13-8  pm  window 
[9,10,111.  The  nature  of  the  continuum  based  on  these  measurements 
was  uncertain.  It  could  be  due  to  far  wings  of  strong  absorption 
bands  or  scattering  and  absorption  by  part icul ates. 

In  an  effort  to  determine  the  cause  of  continuum  absorption  in 
the  13  to  8  pm  window  Bignell  [12]  in  1963  examined  solar  spectra  while 
testing  the  atmosphere  for  aerosol  amounts  and  studied  C0~  far 
wing  contributions.  He  concluded  that  the  amount  of  continuum 
absorption  observed  could  not  be  accounted  for  by  aerosol  attenu¬ 
ation  or  far  wings  of  CCL.  An  attempt  was  then  made  to  model 
the  continuum  by  far  wings  of  the  bordering  H?0  bands.  The  important 
contribution  from  this  initial  work  was  the  realization  of  the 
major  water  vapor  contributions  to  the  continuum.  A  second  paper 
by  Bignell  [131  in  1970  described  a  careful  examination  of  water 
vapor  absorption  in  the  window  regions  employing  a  multiple  traversal 
absorption  cell  and  grating  spectrometer .  Two  important  facts 
about  the  13-8  window  were  observed  namely  a  large  ratio 
of  water  vapor  self-to-foreign-gas  broadening  ability  and  a  strong 
negative  temperature  dependence.  Both  of  these  results  were  not 
anticipated  based  on  the  far  wing  approaches  of  Bignell' s  1963 
paper.  Also  reported  by  Bignell  was  a  similar  but  much  weaker 
continuum  absorption  in  the  4  pm  region. 

Further  laboratory  measurements  performed  by  Burch  [14,151 
using  a  long  path  absorption  cell  and  grating  spectrometer 
examined  the  10  pm  and  4  pin  regions  in  1970  and  1971,  respectively. 
Confirmation  of  the  large  self-to-foreign  broadening  dependence 
and  negative  temperature  dependence  were  made  in  the  10  p(n  window 
region.  The  4  v<n  continuum  was  more  completely  characterized 
and  verified. 
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In  recent  years  laser  measurements  have  contributed  a  great 
deal  to  the  understanding  of  continuum  absorption.  The  stable 
high  intensity  monochromatic  source  provided  by  the  laser  alleviates 
some  of  the  problems  encountered  in  obtaining  accurate  long  path 
absorption  measurements  which  are  characteristic  of  continuum 
(i.e.,  window)  investigations.  Initial  C0?  laser  measurements 
by  McCoy  (161  in  1969  at  10.6  un  indicated  the  ratio  of  self- 
to-foreign  broadening  ability  to  be  200.  The  value  inferred  from 
line  measurements  is  5  [17].  This  result  confirmed  the  preliminary 
observations  of  Bignell’s  1963  paper  concerning  the  dominance  of 
H.,0  contributions  to  continuum  absorption  and  provided  a  quantitative 
number  for  the  self-to-foreign  broadening  ratio.  Recent  compre¬ 
hensive  studies  at  The  Ohio  State  University  [18,19,201  have  extended 
the  frequency  range  by  examining  many  C0~  laser  lines  using  a 
stablized  laser  with  both  a  multiple  traversal  cell  set  to  1.5 
kilometer  path  length  and  also  with  an  optoacoustic  spectrophone. 

The  results  further  verify  the  general  trend  of  large  self  to 
foreign  broadening  ratios  and  negative  temperature  dependence. 

I 

Arefev  and  Dianov-Klokov  (21]  have  also  examined  the  tempera¬ 
ture  dependence  In  the  10.6  pm  region  using  a  C0p  laser  and  long 
path  cell.  A  negative  1.7%  change  in  absorption  per  degree  Kelvin 
was  observed,  over  a  284  to  353  K  range. 

From  extended  temperature  measurements  to  473  K  using  a 
laser  diode  operating  at  1203  wavenumbers  Montgomery  [22]  has 
observed  that  above  398  K  a  leveling  of  the  negative  temperature 
dependence  occurs.  To  explain  the  result  a  competitive  effect 
between  dimer  and  far  wing  absorption  was  suggested.  According 
to  Montgomery,  temperatures  below  398  K  dimer  absorption  dominates 
and  above  398  K  far  wing  absorption  dominates.  Because  these 
measurements  were  made  around  8  pm,  close  to  the  band  of  FLO, 
it  is  not  known  if  a  similar  temperature  dependence  will  occur 
at  10  pm. 

Other  popular  lasers  used  in  probing  the  continuum  have 
been  CO  and  DF  which  operate  at  the  5  and  4  pm  regions  respectively. 
Absorption  at  5  pm  features  strong  local  line  absorption  along 
with  continuum  contributions  [23].  Thus  many  of  the  special  con¬ 
tinuum  characteristics  are  masked  by  comparable  local  line  absorption. 
Recent  measurements  at  White  Sands  [24]  with  a  DF  laser  and  1.5  km 
path  lengths  show  continuum  absorption  properties  at  4  pm  similar 
to  those  in  the  10  pm  region.  An  average  self-to-foreign  broadening 
ratio  of  around  90  has  been  measured.  Thus  it  is  possible  that 
the  4  pm  window  continuum  may  be  generated  by  the  same  phenomena 
as  the  10  pm  window. 


Modeling  efforts  have  not  been  totally  successful  in  explain¬ 
ing  the  experimental  observations.  Attempts  by  Burch  [141  to 
model  measured  absorption  in  the  8  -  12  urn  region  using  the  Lorentz 
line  shape  failed  to  predict  the  observed  absorption  magnitude, 
water  vapor  partial  pressure  dependence  or  temperature  dependence. 
This  computation  used  all  absorption  lines  in  the  rotational  and 
^2  bands  of  H~0  (i.e.,  no  bound).  Even  though  the  validity  of 
tne  Lorentz  line  shape  does  not  extend  beyond  5-10  cm"  ,  the 
result  has  been  a  de-emphasis  of  far  wing  models  in  general. 

This  has  prompted  the  search  for  other  mechanisms  of  continuum 
absorption.  However,  sophisticated  far  wing  theories  are  being 
applied  to  the  HpO  continuum  by  a  few  researchers. 

A  major  candidate  to  describe  continuum  absorption  has  been 
the  water  vapor  dimer  [25,261.  The  strength  of  this  approach 
is  the  strong  negative  temperature  dependence  exhibited  by  dimer 
absorption.  A  similar  candidate  is  the  water  vapor  cluster  [1,271 
which  extends  the  dimer  concept  to  larger  formations.  Because 
monomer,  dimer  and  cluster  formations  are  present  in  any  sample 
of  water  vapor  it  is  difficult  to  determine  the  most  important 
continuum  mechanism  involved. 

Advanced  far  wing  models  are  being  used  by  Fomin  and  Zuev 
[28,29  |  to  explain  experimental  measurements  of  the  HpO  and  COp 
continuums.  Significant  improvements  over  the  Lorentz  line  shape 
far  wing  are  demonstrated  for  the  HpO  partial  pressure  dependence 
and  frequency  dependence.  However  a  lack  of  detailed  information 
on  the  nature  of  this  (i.e.,  Fomin's)  line  shape  formula  make 
it  difficult  to  properly  evaluate  its  performance.  Nonetheless, 
the  success  of  their  far  wing  approach  has  encouraged  our  efforts 
to  develop  a  total  line  shape  formula  for  HpO. 

The  amplitude  of  the  far  wings  of  water  vapor  spectral  lines 
are  difficult  to  determine  both  theoretically  and  experimentally. 
Figures  2  and  3  show  the  frequency  dependence  of  the  10  |jm  and 
4  urn  continuum  regions  respectively.  The  observed  frequency 
dependence  is  consistent  with  a  far  wing  extension  of  the  bordering 
HpO  bands.  The  strong  HpO  rotational  band  should  dominate  in 
tne  10  |,m  region  from  the  low  frequency  side  as  Figure  2  indicates. 
The  4  region  exhibits  more  balance  in  the-jStrength  of  the  border¬ 
ing  bands.  The  Up  fundamental  near  1600  cm’  contributes  from 
the^low  frequency  side  and  the  v,  and  v,  fundamentals  near  3700 
cm_.|  contribute  from  the  high  frequency  s ide.  A  minimum  at  2600 
cm  occurs  in  the  continuum  due  to  the  nearly  equal  strengths 
of  the  Up  and  the  v,  and  bands.  Based  on  the  frequency  behavior 
of  the  HpO  continuum  far  wings  must  contribute  some  portion  of 
the  observed  absorption,  with  perhaps  dimer  and  cluster  absorption 
accounting  for  the  rest. 
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Figure  2.  Water  vapor  absorption  C°(v)  at  T=296K  as  a 
function  of  frequency  v  (wavelength  X) 
in  the  8-12  pm  region. 


Far  wings  of  water  vapor  absorption  lines  are  also  significant 
in  the  fundamental  bands.  The  FLO  spectrum  in  Figure  1  shows 
wing  absorption  underneath  spectral  structure  at  the  band  edges 
of  the  v2  band  at  5  pm  and  7  to  8  pm.  Far  wings  of  absorption 
lines  within  the  band  can  be  studied  in  such  regions.  Once  the 
nature  of  the  far  wings  is  determined  in  the  near  band  region 
an  extrapolation  in  frequency  away  from  the  band  and  into  the 
windows  can  be  performed.  Thus  the  characteristics  of  far  wing 
absorption  can  be  determined  in  the  window  regions.  To  do  this 
a  proper  model  for  far  wing  absorption  must  be  developed  and  in¬ 
corporated  into  a  total  line  shape  theory,  from  line  center  to 
far  wing  of  a  single  absorption  line.  The  total  line  shape  must 
be  normalized,  thus  insuring  conservation  of  energy  and  proper 
magnitude  of  absorption.  To  achieve  this  goal  an  in  depth  analysis 
of  the  interactions  of  photons  and  gases  is  required. 


FREQUENCY  (cm- 


750  torr  of  air. 


The  purpose  of  this  study  is  to  demonstrate  the  importance 
of  far  wing  phenomena  in  characterizing  H?0  continuum  absorption. 
Chapter  2  attempts  to  develop  a  valid  total  line  shape  for  water 
vapor-ni trogen  interactions  under  tropospheric  conditions.  Chapter 
III  describes  the  experimental  apparatus  and  procedures  used  in 
making  continuum  measurements  on  White-type  absorption  cells. 
Chapter  IV  presents  the  experimental  results  of  this  study  and 
a  comprehensive  review  of  other  experimental  efforts.  Using  this 
data  base  parameters  of  the  far  wing  component  of  the  total  line 
shape  are  determined  from  near  band  experimental  data.  Then 
a  comparison  of  experimental  and  theoretical  results  is  made 
displaying  remarkable  agreement.  The  total  line  shape  calculations 
are  applied  to  the  4  pm,  5  vim,  10  p,  17  pm,  and  millimeter  windows 
under  tropospheric  conditions. 

Appendix  B  contains,water  vapor  spectra  inthe  window  regions 
from  850  cmT  to  1282  cm’  ,  1800  cm’  to.2054  cm,  and  2500  cm’1 
to  2862  cm’  .  The  resolution  is  .05  cm’  . 
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CHAPTER  II 

THEORY  OF  SPECTRAL  LINES 


A.  Introduction 

A  general  expression  for  the  absorption  coefficient  is 
developed  in  this  chapter.  Section  B  derives  a  general  expression 
from  first  principles  containing  line  strength  and  line  shape 
information.  Section  C  develops  the  line  strength  resulting  in 
a  standard  formula  with  some  correction  terms.  Section  D  separately 
develops  line  shape  expressions  valid  near  line  center  and  the 
far  wing  under  the  adiabatic  assumption  and  binary  collision  approxi¬ 
mation.  Section  E  presents  a  technique  for  constructing  the  total 
line  shape  for  H?0-N?  gaseous  mixtures  and  a  careful  analysis 
of  the  intermolecular  potential  functions  for  H^O-H^O  and  H^O- 
N„  collisions.  Knowledge  of  the  intermolecular  potentials  Is 
crucial  in  predicting  the  far  wing  behavior  of  a  line  shape. 

The  important  equations  of  this  chapter  will  be  enclosed  in 
rectangles. 

In  order  to  examine  the  attenuation  of  monochromatic  radiation 
per  unit  path  length,  consider  a  medium  (Figure  4)  on  which  plane 
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wave  electromagnetic  radiation  is  incident.  Let  "I"  be  the  intensity 
at  any  point  along  the  path  and  "H"  the  path  length.  The  medium 
is  divided  into  a  number  of  segments  of  length  d  perpendicular 
to  the  direction  of  the  radiation.  The  change  in  intensity,  dl, 
per  unit  length,  dH,  is  proportional  to  some  constant  times  the 
intensity  plus  another  constant  independent  of  intensity  (for 
a  homogeneous  medium).  This  may  be  expressed  as. 


The  parameter  o  is  the  extinction  coefficient  and  e  is  the  emission 
coefficient.  In  the  following  analysis  the  emission  coefficient 
is  not  considered.  The  emission  coefficient  represents  blackbody 
radiation  of  the  medium  and  future  arguments  will  justify  its 
elimination.  The  physical  processes  which  contribute  to  the  ex¬ 
tinction  coefficient  can  be  divided  into  two  general  classes, 
absorption  and  scattering.  Thus, 

a  =  k  +  s  (2) 

where  k  =  absorption  coefficient 
s  =  scattering  coefficient. 


The  dominant  coefficient  will  be  determined  by  the  nature  of  the 
medium  and  by  the  radiation  frequency.  This  study  is  confined 
to  gaseous  media  and  the  infrared  spectrum.  This  condition  basically 
eliminates  scattering  (because  a«X),  which  in  the  infrared  becomes 
important  for  aerosols,  clouds,  fog,  and  various  other  particulates. 
As  the  next  section  will  show,  scattering  is  at  least  a  two  photon 
process  where  absorption  can  be  a  single  photon  interaction. 

In  a  perturbation  expansion  the  single  photon  process  will  be 
favored. 

Lambert's  law  for  absorption  of  radiation  is  obtained. 


(3) 


The  solution  of  this  simple  differential  equation  leads  to  the 
definition  of  another  important  parameter,  T,  the  transmittance. 


T  = 


where 


(4) 


The  absorption  coefficient,  k,  represents  the  molecular 
interaction  with  an  electromagnetic  field  and  other  molecular 
systems,  resulting  in  the  absorption  of  a  photon  and  the  excitation 
of  the  absorbing  molecular  system.  The  analysis  of  such  phen¬ 
omena  requires  an  understanding  of  molecular  spectroscopy,  quantum 
theory  of  radiation  and  intermolecul ar  forces.  For  only  one  absorbing 
molecule  and  one  perturbing  foreign  molecule  the  absorption  coef¬ 
ficient  may  be  functionally  expressed  as 


k(Ma,Mf,v)  •  Stu(Ma,v0)j(,-v0.Ma,Mf;a(Ma,Nf,v0)) 


where 


j(Av;ot) 


-  represents  the  absorbing  molecule 

-  represents  the  foreign  molecule 

-  the  linecenter  wavenumber 

-  the  line  strength  for  transition  between 
it  and  u  states 

-  the  1 ine  shape 

-  the  half  width  at  half  maximum 


Note  that  the  line  strength  is  a  function  of  the  absorbing 
molecule  only  (to  first  order).  It  is  convenient  to  normalize 
j(Av;«)  by  requiring. 


M«a,Hf.v)dv  -  Slu(«a) 


ao 

j(v-vo’Ma’Mf;a(Ma’Mf’vo))dv  =  1 

/ 


The  parameters  M  and  M,  are  themselves  multi-variable  functions 
of  quantities  affecting' the  properties  of  the  associated  molecular 
systems.  They  are  used  here  to  emphasize  this  fact  and,  having 
served  their  purpose,  shall  be  eliminated  from  the  notation  and 
replaced  by  more  explicit  variables. 


The  description  of  the  absorption  coefficient,  k(v),  can  now 
be  broken  down  into  two  parts,  line  strength  and  line  shape  function. 
These  two  quantities  are  natural  choices  for  experimentally  charac¬ 
terizing  the  absorption  of  a  single  line.  Also  the  absorption  line 
parameter  compilation  of  AFGL  [5]  uses  this  formulation.  Therefore 
the  goal  of  the  following  theoretical  development  is  to  find  standard 
mathematical  forms  of  lint*  strength  and  line  shape. 

L> .  ,  im  Electrodynamics  of  Gases  and  Photons 

In  the  quantum  mechanical  model  representing  gases  and  photons 
the  Hamiltonian  is  complex;  because  of  this,  a  careful  description 
is  necessary.  The  total  Hamiltonian  may  be  written 


"Tod1* 


'a,  *  Hb  *  Hp  *  Ha,p(t)  *  Hc(t> 


where 


H 

P 

H  (t) 
»c(t) 


-  Hamiltonian  for  the  absorbing  molecule 

-  Hamiltonian  for  the  broadening  molecule 

-  Hamiltonian  for  the  photon  field 

Interaction  Hamiltonian  between  absorbing 
molecule  and  photon  field 

-  Collision  Hamiltonian  between  absorbing 
molecule  and  broadening  molecules. 


Time  dependent  terms  are  explicitly  indicated.  Otherwise  the 
separate  Hamiltonians  are  time  independent.  The  first  three  terms 
in  the  total  Hamiltonian  form  the  unperturbed  or  zero  order  Hamil¬ 
tonian,  H  .  The  last  two  form  the  interaction  Hamiltonian.  A 
description  of  each  term  in  the  total  Hamiltonian  follows. 

H^  represents  the  absorbing  molecule.  In  a  gas  sample 

many  molecules  will  be  of  the  same  type;  the  subscript  1  indicates 
the  absorbing  molecules  while  all  others  are  perturbers.  The 
general  unperturbed  Hamiltonian  for  absorbing  molecules  is 


Ha  =  HaR0T  +  HaVIB  +  HaEL  +  HaTRN 


where 


-  describes  the  rotational  energy  levels 

-  describes  the  vibrational  energy  levels 

-  describes  the  electronic  energy  levels 

-  describes  the  translational  energy  of  the 
total  molecule. 


1? 


H  ...  will  not  be  directly  included  in  the  analysis,  but  will 
be  accounted  for  by  properly  averaging  over  velocities  of  the 
molecules. 

H,  is  a  many-particle  Hamiltonian  representing  all  perturbing 
molecules.  This  includes  foreign  molecules  and  molecules  of  the 
same  type  as  the  absorber.  It  may  be  written 


■  l 


K  "k 

y  y  h. 

k=l  m=l  bkm 


N  is  the  number  of  absorbing  molecules,  K  is  the  number  of  types 
of  foreign  molecules,  and  N,  is  the  number  of  foreign  molecules 
of  the  k"  type.  Each  Hamiltonian  is  structured  in  the  same  way 
as  Equation  (8);  that  is  the  rotational,  vibrational,  electronic, 
and  translat ional  energies  are  represented. 

H  is  the  free  space  photon  Hamiltonian.  Many  excellent 
texts  discuss  this  term  ['ll, .871.  Simply  stated,  it  is 

II  £  flu,  (q)  (affq)  a  (q)  +  l) .  (10) 

P  qo  0  11  ? 

a^fq)  is  the  photon  creation  operator  and  aQ(q)  is  the  photon 
annihilation  operator  for  a  particular  polarisation  o  and  wavenumber 
q  corresponding  to  the  angular  frequency  w=h|q| . 

H  (t)  is  the  interaction  Hamiltonian  between  the  absorbing 
a  jp 

molecule  and  photon  field.  It  can  be  shown  that  [871 


Ha,p(t)  °*{r’t>  °  p) 


♦-S*(F,t)  *S(r,t)>  (u(t)  -  u(T))  (1 

?mc/  J 

p  is  the  particle  momentum  operator.  fi(r,t)  is  the  photon  field 
magnetic  vector  potential  operator. 

>/2  f  _ i  (q  °  r-utl 


^(r.t)  =  \a0(d)  e0(q)  e 


+  a.t(q)  e  (q)  e 


•i(q  °  r-mtl 


o  a 


13 


e  is  an  orthonormal  vector  perpendicular  to  the  free  space  photon 
propagation  representing  the  polarization.  T  is  the  observation 
time  or  the  time  period  during  which  the  molecules  and  photons 
can  interact. 


Hc(t)  models  the  intermolecul ar  potentials  due  to  neighboring 
molecules  surrounding  the  absorbing  molecule.  It  is  the  least 
understood  term  in  the  Hamiltonian  and  yet  when  considering  line 
shapes  it  is  the  most  important.  A  general  many-body  expression 
can  be  written  as 


Hc(t) 


t 

3=2 


Ha  a  (t) 

Vl 


k  Nk 

y  y  Hh  (t) 

k=l  m=2  bmkal 


(13) 


Each  term  is  a  binary  interaction  of  the  external  fields  of  the 
molecule.  Figure  5  illustrates  such  an  interaction.  For  mutual 
neutral  molecules  with  permanent  electrostatic  potentials  a  single 
Hamiltonian  in  Hc(t)  may  be  written  [33] 

Hbmkcii^R^t^’rmk’rl^  ~  ^Electrostatic  +  VMagnetostatic />Rot 


-  V 


Dispersive 


-  V 


Inductive 


+  V 


Repulsive’ 


(14) 


R(t)  i  ,  the  intermo2ecular  distance  between  two  moving  molecules. 

The  quantities  rmk , l  locate  the  charged  particles  in  the  internal 
molecular  coordinate  system  for  the  broadener  and  the  absorber, 
respectively.  The  magnetostatic  potential  function  is  significantly 
smaller  than  the  electrostatic  potential  so  it  is  ignored.  A  brief 
description  of  the  potential  terms  will  be  given. 

A  multipole  expansion  of  the  electrostatic  potential  is 

vPP<Fmk-Fl>|  VDQ< Snk  -71  >  W Fmk  -F1  )*vQQ<Fmk  -F1  > 
Electrostatic’  R3(t)  R4(t)  „5(t) 


+ 


(15) 


where 
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vnn^'mk,ri) 

Wjnk’^j) 

Wjnk’^ 

VQQ(rmk,r^ 


Dipole-dipole  interaction  term 
Dipole-quadrupole  interaction  term 
Dipole-octupole  interaction  term 
Quadrupole-quadrupole  interaction  term 


and 


VDD(rmk’rl) 

VDQ^rmk,rl^ 

VD0(?mk*?l) 

'Wrmk,rV 


Mmklll  f^mk’Ql^ 

Mmk  ^1  9l(Qmk’V  +  ^mk  M1  92(Qmk,nl) 
ymk  *1  hl^mk,JV  +  *mk  M1  h2^mk,SV 
%ik  ^1  ^nmk*nl^ 


Pm^  and  p i  are  the  dipole  moments  of  broadener  and  absorber,  re¬ 
spectively.  IJ  is  the  mean  quadrupole  moment  and  R  is  the  mean 
octupole  moment.  The  functions  of  ft  ,  and  ft,  represent  the  angular 
dependence  of  the  potential  f unct ions  and  whether  the  potential 
function  is  attractive  or  repulsive  at  the  generalized  orientation 
angles  ft  .  and  ft^  in  a  space-fixed  coordinate  system.  The  molecules 
are  rotating  during  the  collision  with  the  amount  of  rotation 
observed  during  the  collision  depending  on  the  collision  duration 
time.  The  forces  induced  by  the  electrostatic  potential  will 
always  try  to  orient  the  molecules  in  a  minimum  energy  configuration 
resulting  in  a  net  attractive  force.  An  average  over  this  rotation 
is  performed  to  account  for  the  molecule's  ability  to  achieve  this 
preferred  orientation;  it  is  denoted  by  < 


Rot* 


<VE1>, 


/' 


-Vr 


dft  VE1  e  /  P(n0’td’fmk’f 1^ 


Rot 


f 


(16) 


HO  /VE1/kT  nfO  t  fR  fR\ 

e  P^o^d’^nk’^ 


where  dft  is  the  differential  element  for  all  orientation  angles 
and  p(  » t  ^ ,  f  » f  7 )  's  a  pulse  function  with  unit  amplitude  and 
pulse  duration  depending  on  ftQ,  the  initial  orientation,  t^,  the 
collision  duration  time  and  fR^  and  f*j,  the  rotational  frequencies 
of  the  broadening  molecule  and  absorbing  molecule,  respectively. 


I 
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For  collisions  with  large  impact  parameters  the  collision  duration 
time  is  very  brief,  thus 


C^[))  X  (17) 

L  Rot 

For  collisions  with  small  impact  parameters  the  electrostatic 
potentials  can  be  so  strong  that  it  is  no  longer  a  small  perturbation 
to  the  rotational  motion.  This  happens  in  the  formati  of  liquids 
and  solids.  The  molecule's  ability  to  rotate  is  inhibited  as  it 
takes  a  relatively  fixed  orientation.  Thus  the  average  over  rotations 
breaks  down.  In  fact  the  formulation  of  the  unperturbed  Hamiltonian 
is  no  longer  valid.  This  is  an  important  point  for  water  vapor 
and  will  be  discussed  in  more  detail  later. 


Not  all  molecules  have  permanent  dipole  moments,  and  thus 
other  types  of  long  range  potentials  are  of  interest.  The  most 
important  are  jgduction  and  dispersion,  both  with  leading  terms 
of  the  form  CR~°(t).  Induction  potentials,  ^ j nduCt i ve *  rePresent 

the  interaction  of  a  permanent  electric  multipole  of  one  molecule 
with  an  induced  electric  multipole  of  another  molecule  with  perhaps 
no  permanent  multipole.  The  resulting  potential  is  attractive. 
Dispersion  potentials  arise  from  the  random  motion  (within  the 
bounds  of  the  uncertainty  principle)  of  electrons  about  the  nuclei 
generating  instantaneous  multipoles.  These  brief  multipoles  can 
induce  multipoles  in  other  molecules  resulting  in  an  attractive 
potential.  For  a  detailed  description  of  these  potentials  see 
Margenau  and  Kestner  [33].  Approximate  forms  are  given  below 
[34]. 


Inductive 


Dispersive 


E 


2  -  .  2  - 
P 1  “mk  Pmk 


al_ 


R6(t) 

-  q2  2 
_  fl  e_<r  > 

2  R6( t) 


(18) 

(19) 


where 


a  -  mean  electric  polarizability 
e  -  charge  of  an  electron. 


The  discussion  of  H  (t)  thus  far  has  involved  the  long  range 
potentials.  Intuitively.these  are  the  most  important  potentials 
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to  consider  for  gases.  The  short  range  potentials,  VpepU]sjve» 

are  repulsive  in  nature  such  as  experienced  by  the  electron  clouds 
of  very  close  molecules.  The  combination  of  the  long  range  attractive 
potentials  and  the  short  range  repulsive  potentials  determine 
intermolecular  bonds  which  result  in  possible  polymer  formations. 

In  some  gases  the  dimer  is  by  far  the  most  common  polymer.  Knowledge 
of  bond  lengths  and  strengths  are  helpful  in  determining0t.he  nature 
of  intermolecular  potentials  for  small  separations  ( 2-7  A).  Figure 
6  shows  the  form  of  a  general  intermolecular  potential. 


Figure  6.  The  intermolecular  potential 


The  development  of  spectral  line  intensity  by  Anderson, 

Tsao  and  Curnutte  [35,36]  will  be  followed  to  obtain  a  transition 
rate  of  the  system.  The .probabi 1 ity  that  the  system^will  change 
from  an  initial  state  | s 1 ( t ) >  to  some  final  state  <s'(t)|  is 

P._^f  =  |<sf(t)|  si(t)>|2  .  (20) 

The  zero  order  eigenfunction  of  the  system  is  a  product  of  the 
eigenfunctions  corresponding  to  each  individual  zero  order  Hamil¬ 
tonian  composing  the  total  quantum  system. 
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I 


|s(t)>=|a1(t)a?(t)...aN  (t)b11(t)...bk  (t)...bKN  ( t) . . .n-  ( t) . . .> 

a  k  4 

(21) 


The  notation  follows  from  the  definition  of  the  Hamiltonians. 

The  time  dependence  of  the  initial  and  final  eigenstates  can  be 
characterized  by  the  appropriate  time-development  operator  (scattering 
operator) 


|si(t)>  =  ST()t(t)|  si>  (22) 

and 

<  s f ( t ) |  =  <sf |S”^(t) .  (23) 


(24) 


(25) 


H  (t)  lacks  only  the  interaction  between  field  and  particle  since 
tn£  final  state  comes  after  the  photon  interaction.  Therefore, 

pi.f  •  l«fls«<‘>  (26) 

Let 

*<*>  -  S>>  sIot<‘>  (?7) 


A  convenient  form  for  the  operators  of  concern  is 


TS  HTot(t')dt' 


sTot<‘>  ■  e 


So)(t)  *  e 


t 

is  j  Hoc<«dt' 


where 


"oc*1’  *  Ho  *  Hc<*>  • 


then 


is  *(t)  ■  (s;J(t)  Hyt)  soc(tn  v(t) 

\ 
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(28) 

I 


I 


where  the  relations 


,fl  S>>  *  '  S>>V‘> 

STot^  =  HTot^STot^ 
have  been  used.  Equation  (27)  leads  to 

V(t)  =  e*p  ^  J  S«(f)H4ip(f)S0c(f)c«j  •  (29) 

A  perturbation  expansion  of  the  above  result  gives 

t  f  t(n-T) 

TO  Jdt-j  dt<"»(s;J(f)H  p(t')S0C(f)) 

0  0  0 

‘Soc<t,,>Ha,p't">S0C<t"’>---(Soc(t<n>)Haip(t(">)Soc(t(">)). 

(30) 


Terms  for  n>2  represent  multiple  photon  processes  such  as  Raman 
scattering,  Rayleigh  scattering,  second  harmonic  generation,  etc. 
For  the  Photon  density  levels  and  frequencies  of  interest  these 
higher  order  terms  are  unimportant.  V(t)  to  first  order  is 


V(t)  =  1  + 


t 


1 


dt'  s«<t,>Ha)plt'>S„c<f> 


(31) 


Recall  also  that  H,  (t)  contains  two  terms,  one  involving  a  sinqle 
a  ^  p '  ’  3  3 

photon  interaction  (P  °  A)  and  the  other  a  two  photon  process 
(A  ®  A).  For  the  same  reasons  as  above  only  the  P  °  A  term  needs  to 
be  retained.  Also,  because  the  eigenstate  has  the  particle  and 
photon  states  independent  of  one  another  the  interaction  Hamiltonian 
can  be  written  using  Equations  (11)  and  (12)  as  follows. 


1/2 


V'>  ■  <-  li)  (I2#1)  (*„(5)P  •  e  (q)e1'>°'V,“t 


1/2 


+  (^T}  (a*{q)  P  °  ^(q)e‘iq°r)eiu,t)(u(t)-u(T)) 

(32) 


\pM  = 


Hab  loot  ^  ^em 


a  1 P 


(u(t)  -  u(T) ) . 


(33) 


Equation  (33)  has  two  parts,  one  for  absorption  and  one 
for  emission  of  a  photon.  For  absorption  the  matrix  element  becomes 


<sf|vab(t)|si;>=  ■!*  j  dt'<sf|s^(f)H|bp  soc(f)| 


(34) 


and  for  emission 


<sf|vem(t)|si>=  df<ff|s^{t)H“  Soc(f)  ls4> 


Ju)t' 


(35) 


The  zero  order  part  of  the  time  development  operator  is  removed 
by  observing 


Soc{t)|s1>=  exp^Tfi  Eil)  sc(^|s1> 


(36) 


Csf|S^(t)  =  exp(^  Eft)<'sf|S-1(t) 


(37) 


Also,  using  ti^  =  Ef-E.j  and  Aujfu>-u)o  the  matrix  elements  become 
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<sf|Vab(t)|si>-  ^  dt 1  e 


and 


'iAwt'  \sf  |S^1(t')H^pSc(f)|si> 

(38) 


Vsf|Vem(t)|sl>  =  I  dt'  e 1  Atut  ‘<sf  |  Si 1  ( t ' )  H^m  S  ( t ' )  |  s  ^  •  ( 39 ) 


a-,p 


It  is  only  necessary  to  analyze  Equation  (38)  since  the  emission 
matrix  element  will  develop  analogously.  This  is  a  manifestation 
of  detailed  balancing. 

Substituting  Equation  (38)  into  Equation  (26)  yields 


P4 


j-f  ti2 


dt’  e 


,<5f|S;1(t)H^pSc(t,)|si> 


(40) 


The  general  probability,  P,  is  obtained  by  summing  P-^  over  the 
final  states  and  averaging  over  the  initial  states 


p*<Ip^f>i  • 

(41) 

probability  of  a  transition  becomes 

T  T 

^  Jdt J dt"  e"lAw(t"_t,)Tr 

>Ha^pSc<f>] 

0  0 

(42) 

where  p  is  the  density  matrix  generated  by  the  initial  states. 
Since  tRe  operators  in  a  trace  can  be  cyclically  permuted  the 
trace  in  Equation  (42)  becomes 

Tr[eoHaVS^(t")Sclt'),HV(Sc(t")S'1(t'))  ’ 


It  is  convenient  to  observe  that 
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and 


S‘](f)  Sc(f)  =  S~7(t) 

Sc(t")  S~C](V)  -  Sc(x)  ( 

where  T»t"-t*.  Tsao  and  Curnutte  [36]  maintain  these  relations 
to  be  true  when  averaged  over  all  molecular  collision  paths 
This  simplifies  the  trace  to  the  form 


F(t)  .  rr[p^ps-'(T)„*psc(T)]  .  I 

Note  that  F { t)  =  F  ( -t)  (where  *  denotes  complex  conjugate)  and 
transform  the  variables  of  integration  from  t"  and  t1  to  t  and 
t  .  Figure  7  illustrates  this  change  of  variables.  Therefore, 


T  T-t  0  T 

p  =  j  dT  F(t)  J  dt '  +  J  dt  F(t)  Jdt 


f 

P  »  ^  Rejdr  (T-t)  f(t) 


For  a  cw  probe  photon  field  the  observation  time  is  extremely 
long  compared  to  interaction  times. 

It  is  assumed  therefore  that  T»t  and  P  becomes 


i 

P  -  ^  Re  f  dx  F ( t) 
J 


The  assumption  that  T»t  will  in  practice  become  invalid  for  100 
?1,n-^60nd  or  7*ss  Pulse  durations  of  the  photon  field  (i.e., 
|<I0  seconds).  The  transition  rate,  r,  is  the  probability 
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Figure  7.  Change  of  the  integration  variables. 


of  observing  a  transition  per  unit  of  time  or  r  =  y.  The  transition 

is  between  two  energy  levels  designated  upper  and  lower.  For 
absorption,  the  lower  energy  is  the  initial  state.  Thus  the 
transition  rate  can  be  written  ,  and  Equation  (47)  becomes 


•  (48) 


The  observation  time  has  been  extended  to  infinity  with  very  little 
loss  in  accuracy.  Equation  (48)  forms  the  basis  for  many  modern 
line  shape  theories  [37,38,39,40,41]  which  have  followed  the  work  of 
Anderson,  Tsao,  and  Curnutte. 
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Equation  (48)  is  written  in  a  very  compact  form  and  will 
be  expanded  by  summing  over  the  initial  states  of  the  trace  and 
a  complete  set  of  final  states,  i.e. 


where 


%  -  'V1 

g^  V  pVibpRot 


degeneracy  of  lower  level 

density  matrix  over  vibrational 
and  rotational  states 


density  matrix  for  the  photon 
field 

density  matrix  over  vibrational 
and  rotational  states  of 
broadening  molecule. 


The  "a"  quantum  numbers  represent  all  other  quantum  numbers  necessary 

to  describe  the  absorbing  golecule  besides  the  rotational  quantum 

numbers  J  and  M.  s^  and  s°  represent  all  states  and  necessary 

quantum  numbers  for  the  broadening  molecules.  To  abbreviate  the 

notation  let  fc=a9J9Mff  and  u=a  J  M  and,  real  izing  H,  is  independent 

.  b  u  u  u  u  alp 

of  s9  or  s.  we  have 

**  u 


?Re  (I  T.  Y. 

alJ*.  "Vmu  qt°ll 


q<y 


I"  "qau>  b  s 


q  ,cr 
Mu’  u 


LKbwbb 
b  b  si  sesu 

Vi 


j  d,  «-,4>  ^  Vu1  Sc<T)llst  V>  ' 

0 
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r£u  =  ^  Re 


Reducing  the  sums  for  nontrivial  matrix  elements 
2* 


Z  H  ^  Z  Z  O  J  H V1  “n%u>  Z  Ps 

N  1  *  M  M  -  •  - 


a£’J£  M£’Mu  Vfc 
V°u 


» f  *  «,  \js;'(t)4  sc(t]us\o1>! 


(51) 


The  subscripts  "u"  and  "£"  have  been  dropped  from  sb  and  sb  in 

the  above  formula  since  the  Kronecker-6  requires  £=u.  By  inserting 

unit  operators  between  the.  time  development  operators  and  the 

interaction  Hamiltonian  H,  the  following  relation  is  obtained 

a^p  3 


<“S\0U  Is''  Sc<T)llsV>" 

VUlSc’(T)lu'sb' 

x4,sl>“  n^lSc<^l‘st>V>  (52) 

Since  S  is  independent  of  the  photon  field  and  Hab  1S  independent 

c  b'  b"  alq 

of  the  broadening  molecules  s  =s  ,  n'  =  n„_ ,  and  n'  „  =  n„ 

*  qo£  qo£’  qau  qau 

Further  the  adiabatic  assumption  [36,37]  is  made  requiring  u=u', 
b  b* 

s  =s  ,  and  £=£'.  Therefore,  Equation  (52)  becomes 

<usb|S-'(T)|usbXunq(jn|H^p|)lnqij)l)<tsb|Sc(T)|tsb>.  (53) 


Z. .  «*' 


u'sb' 


n'  £' 
qou 

sb"n' 

qa£ 
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I 


Substituting  this  result  into  Equation  (51)  we  obtain. 


1  r 


1 


i 


L 


r«u  =  17  Z  L  ppKuVulHfnltn  Vl? 


f)  —  -  ^0 

a^,  *  MtMu  q^ 

V°u 
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x  ]T  p  b  <usb|S-1(x)  |usbXesb|Sc(T)|!tsb>  .  (54) 


Now  examine  a  single  matrix  element  of  the  time  development  operator 
in  terms  of  each  collision  Hamiltonian  between  the  absorber  and 
all  broadeners 


Na  "l  Nk 


itsb|s  (t) Ins'S  - <tsb I TT  s  n  s  ...rf  s  |esbS 
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(55) 


Assuming  binary  collisions, 

slsb|Sc(T)|tsb>  -<*sa|ScaiaHsa>  a  Vlsb,|Sca]bi|tsb|> 


bk  bk  ^ 

•sts  lSca,bklts  > 


(56) 


Scaja’  ^ca^bp  ^ca^  rePresent  quantities  averaged  over 

velocity  and  impact  parameter,  and  summed  over  collisional  events. 
Analogously,  it  can  be  stated. 
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(57) 
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The  sum  over  broadening  states  in  Equation  (54)  can  be  expanded 
to  be 


The  symbol  <  >R  indicates  the  average  over  impact  parameter 

o’  ’  o 

and  velocity  and  sum  over  collision  events.  It  is  written  mathe¬ 
matically,  for  arbitrary  A, 


<A> 


-j  -  ( MV  /k T )  r  5- — 

vdve  I  RodRoLA 

RMin  to 


where 


(59) 


V  -  vo 1 ume 

p  -  reduced  mass  of  binary  system 

T  -  temperature 

k  -  Boltzmann's  constant. 


Also  because  N  »1  then  N  -1~N  this  fact  is  used  in  Equation  (58). 
It  is  convenient  to  define 
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The  transition  rate  takes  on  the  familiar  form  of  Fermi's  Golden 
Rule. 
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(62) 


The  sums  over  a»  and  J.  simply  add  the  effects  of  all  possible 
transitions.  Tne  terni  in  parenthesis  in  Equation  (62)  leads  to 
the  absorption  contribution  to  the  line  strength  for  a  single 
transition.  To  complete  the  line  strength  emission  processes 
must  be  included.  The  transition  rate  for  emission  is  (Hc(t)=0) 
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(63) 


The  frequency  dependence  of  the  transition  rate  is  contained  in 

the  integral  expression  and  this  becomes  the  line  shape  for  collision- 

broadened  lines.  Thus 

CD 

j^(u))  =  ^  Re  j  dt  e‘iAu)T  C(t)  .  (64) 

o 

C(t)  is  interpreted  as  an  autocorrelation  function  between  the 
state  of  the  atom  at  x=0  and  at  time  t.  It  is  symmetric  about 
t=0,  monotonical ly  decaying  as  time  increases.  Using  the  symmetric 
property  of  C(t),  Equation  (64),  can  be  re-expressed  as 
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R  U0oT 

Thus,  j  (w)  and  e  C(-r).are  Fourier  transform  pairs.  To  generalize 

the  concept  further,  let  j 1 (to )  be  the  Hilbert  transform  companion 
of  j*(u)).  Define  c 


jcM  =  r 


C(t)>  =  j,(u>)  +  ij*(u>)  . 


j 1  (to )  represents  refractive  index  effects  of  the  collision  process. 
Because  the  interest  of  this  study  is  centered  upon  line  shapes 
(refractive  index  changes  are  small  for  gases)  the  superscript 
R  will  be  dropped  and  the  line  shape  function  will  be  denoted 

by 


jcM  =  Re  F' 


no  t 

e  0  C  (  t  ) 


The  line  strength  and  line  shape  equations  will  be  developed 
separately  in  the  following  sections. 


The  Line  Strength 


The  line  strength  of  a  particular  transition  (fixed  a-,J.)  will 
now  be  considered.  For  a  medium  in  thermal  equilibrium  the  rate 
of  transitions  from  particular  upper  to  lower  energy  states  and 
lower  to  upper  energy  states  must  be  equal,  i.e. 


dN  \ 
_Jo 
dt  | 


absorbed 


emitted 


N  is  the  photon  occupation  number,  that  is,  the  number  of  photons 
wrlh  wavenumber  q  and  polarization  a ■  The  quantities  in  Equation 
(68)  are  related  to  the  transition  rate  by 
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where 


(70) 


n  is  the  total  number  of  absorbing  molecules  per  volume.  Using 
Equations  (62)  and  (63)  for  fixed  a  ,J  the  above  relations  are 
written, 
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(72) 

n  and  nau  are  number  absorbing  molecules  in  the  lower 
and  upper  energy  level,  respectively.  To  obtain  the  standard 
form  for  the  line  strength,  the  matrix  element  in  Equation  (71) 
and  Equation  (72)  is  assumed  to  be  independent  of  the  sum  over 
initial  photon  states.  This  is  usually  a  good  approximation  but 
breaks  down  for  pure  rotational  transitions  and  a  correction  factor 
will  be  introduced  to  handle  this  situation.  Since 
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Equations  (71)  and  (72)  become 
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The  above  expressions  are  reduced  using 
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The  sum  over  the  final  photon  states  is  approximated  by 

Z  Z  ”  3  f  q2  dq  dV 
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The  primed  transition  rates  become 
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(81) 
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The  sums  over  the  final  photon  states  effectively  result 
in  an  average  over  the  random  orientations  of  the  molecules  and 
are  only  valid  for  random  orientations  between  the  photon  field 
and  dipoles.  Substituting  Equations  (81)  and  (82)  into  (75)  and 
(76)  we  obtain 
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is  the  Einstein  A  coefficient. 
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Using  Equation  (68)  the  thermal  equilibrium  relation  is  obtained 
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The  Einstein  B  coefficients  are  obtained  by  defining 

2 

L  w  ■ 


v 

ut  3f.2  % 


B  .  =  A 


=  2jl_  1 


M  A 

u’  Jt 


B  =  —  B 
Jtu  g  uit 


and 


u(uo) 


hu)“ 


W- 


2  3  qail 

71  C 


where  cu(co  )  is  the  intensity  of  incident  radiation.  Equation 
(85)  is  now  written  in  the  familiar  form  of  the  Einstein  relation 
[42) 
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Consider  an  incident  radiation  field  of  frequency  gj  upon 
this  medium  in  thermal  equilibrium.  The  radiation  field  will 
upset  the  balance  between  absorption  and  emission,  giving  rise 
to  induced  transitions,  ntr,  in  unit  time  per  unit  volume. 


tr 
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ae 
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Because  spontaneous  emission  is  essentially  isotropic  and  induced 
emission  and  absorption  are  in  the  same  direction  as  the  incident 
radiation,  the  spontaneous  emission  contribution  in  the  direction 
of  incident  radiation  is  small  and  will  be  ignored.  Therefore, 
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Substituting  for  the  Einstein  B  coefficients 
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Using  a  Maxwel 1-Boltzmann  population  distribution,  ntr  becomes 
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Electronic  partition  function 
Vibrational  partition  function 
Rotational  partition  function 
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Recalling  Lambert's  law 


dl 
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kl, 


and  multiplying  both  sides  by  a  differential  of  frequency  yields 

-  dco  =  kldco  .  (91) 


Recall  that  the  spectral  intensity  I  =  cu(w)  represents  the  incident 
energy  per  area  and  also  that  the  negative  change  in  intensity  per 
length  over  an  incremental  band  of  frequencies  is  equal  to  the  net 
number  of  transitions  induced  by  the  radiation  field  in  unit  time 
per  unit  volume  times  the  energy  contained  in  that  incremental 
band  of  frequencies.  Thus, 
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If  k(oj)  is  nonzero  over  an  incremental  frequency  range  where 
u(w)  is  slowly  varying  (blackbody  approximation)  or  if  u(w) 
represents  a  monochromatic  source  (laser  approximation)  where 
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S.  is  the  line  strength  of  a  particular  transition  as  defined  by 
Equation  (6).  The  above  equation  is  valid  for  broadband  or  laser 
sources.  Substituting  Equation  (90)  for  n.  ,  the  desired  expression 
for  the  line  strength  is  obtained,  r 
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from  the  ideal  gas  law,  the  number  of  absorbing  molecules  is 


’a  '  FT 


Using  this  fact  and  converting  from  radian  frequency  to  wavenumbers, 
v,  the  line  strength  becomes 
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The  line  strength  at  a  standard  temperature  and  pressure  is 
denoted  as 
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At  any  temperature  and  pressure. 
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for  the  AFGL  listing  T  =?96k  and  p  =1  atm.  Also  the  basic  definition 
of  the  absorption  coefficient  is  different  for  the  AFGL  listing 
from  that  defined  in  Equation  (3).  To  convert  between  the  AFGL 
line  strength  (denoted  as  )  and  the  line  strength  defined  by 
Equation  (96),  use 


(99) 


S1U  °-n3m  «  «>:V  <7S0>  t  s?u  • 

for  mil)  to  near  infrared  (v  >  1000  cm"')  near- 1 ine-center  calculations 
the  line  strength  defined  by  Equation  (96)  works  very  well.  However 
for  pure  rotational  transitions  corrections  to  the  line  strength 
must  be  made. 

The  first  correction  for  rotational  transitions  originates 
from  the  molecule's  ability  to  align  its  dipole  to  the  electric 
vector  of  the  photon  field.  If  the  sum  over  initial  photon  states 
had  been  kept  at  the  beginning  of  the  development  in  this  section 
the  correction  factor  to  the  line  strength  would  have  been  included. 
This  correction  is  basically  the  Van  V leck-Wei sskopf  line  shape 
|11,44|  which  includes  a  u-  factor  multiplying  the  Lorentz  line 

shape.  This  correction  oniy  applies  to  near  line  center  of  a 
transition  [4.1I  since  the  impact  approximation  was  made. 

The  second  correction  pertains  to  the  Maxwel 1-Boltzmann 
population  distribution  where  the  energy  levels  are  not  fixed. 

This  accounts  for  population  changes  in  the  shifted  energy  levels 
and  how  they  will  affect  the  population  difference  for  that 
transition.  Let  aE.  and  AE  represent  these  level  shifts  of  the 
lower  and  upper  energy  states  respectively.  Then  using 
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gives  a  correction  factor  of 
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The  value  of  AE^  is  in  practice  unknown.  Near  lint'  center  modeling 
attempts  by  Farmer  1451  to  predict  water  vapor  rotational  band 
absorption  ignores  this  term  (AE,j<<kT). 


The  resulting  correction  to  the  standard  line  strength  equation 
becomes 
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The  line  strength  for  far- infrared  lines  (\><1000  cm  )  becomes 


near  1 Ine  center 


(102) 


E  t^lTTT‘1 

..  ,v.  °  (I..  ° «i  -aE./kt  ,,-.-hwkT) 


“kOJ.  e 


d-e  0  ) 


far  wing. 


figures  8  and  9  show  the  population  difference  factors 
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gure  9. 


become  unity  for  v0  ■»  1000  wavenumbers.  Figure  8  illustrates 
the  negative  temperature  dependence  of  fj(v0,T)  for  different 
values  of  \>Q.  Figure  9  demonstrates  the  importance  of  f^Vg.v) 
in  the  blue  wing  (fp>l,  v-v  =Av>0)  and  red  wing  (f?<l ,Av<0)  °for 
different  values  of  vQ. 

To  summarize  the  results  of  this  section  the  line  strength 
wijj  be  presented  in  four  limits,  greater  than  and  less  than  1000 
cm*  ,  near  line  center  and  far  wing.  Notice  that  some  of  the 
corrections  to  the  line  strength  are  frequency  dependent  and  by 
definition  belong  to  the  line  shape  function. 
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D. 


The  L ine  Shape 


The  most  general  classification  of  line  shapes  is  that  of 
homogeneous  and  inhomogeneous  broadening.  Homogeneous  broadening 
means  that  all  molecules  have  the  same  basic  line  shape  character¬ 
istics.  That  is  if  a  line  shape  is  observed  for  a  collection 
of  molecules  the  same  line  shape  will  be  observed  for  each  molecule. 
Examples  are  natural  broadening  (radiation  damping)  and  collision 
broadening.  Inhomogeneous ly  broadened  lines  represent  a  collection 
of  frequency  shifted  homogeneously  broadened  lines.  Thus  each 
molecule's  line  shape  may  be  completely  different  from  the  total 
line  shape  of  a  collection  of  molecules.  Examples  are  Doppler 
broadening,  non-uniform  electric  and  magnetic  fields  in  Stark 
and  Zeeman  effects,  inhomogeneities  in  a  medium  and  crystalline 
strains  and  defects. 

1 .  Homogeneous  line  shape 

The  Fourier  transform  relationship  of  Equation  (67)  typifies 
homogeneous  1  in  -  shapes.  This  is  easily  demonstrated  using  the 
time-bandwidth  product  of  Fourier  transform  theory. 

At  Aw  =  1 . 

Thus  the  width  of  a  line  is  the  same  for  all  molecules. 

The  line  shape  and  correlation  function  are  divided  into 
three  regions  pertaining  to  near-line,  intermediate,  and  far-wing. 
Separate  conditions  determined  by  the  collisional  dynamics  charac¬ 
terize  each  region.  Figure  10  illustrates  the  effects  of  collisions 
at  different  ranges  between  the  broadening  and  absorbing  molecule. 
The  lines  labeled  1  and  2  represent  possible  collisional  paths. 

Table  1  lists  the  characteristics  of  collisions  in  each  region. 

These  conditions  will  be  applied  to  the  correlation  function  to 
extract  the  form  of  the  correlation  function  to  be  applied  to 
the  near  line  center  and  far  wing  regions.  Figure  11  shows  a 
general  correlation  function  and  where  the  different  types  of 
broadening  will  dominate.  The  functional  form  of  the  line  shape 
is  correspondingly  broken  down  into  three  parts. 
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Figure  10.  The  correspondence  between  impact 
parameter  and  line  shape  region. 


Table  1 

Characteristics  of  Molecular  Collisions 


R  >  R  Very  little  effect  on  absorbing  molecule 

d 

R  >  R  >  Rb  Interruption  broadening  1461  -  near  line 

a  center  phenomena 

1.  Brief  duration  of  collision 

2.  Time  between  collisions  is  long 
compared  to  duration  of  collisions 

3.  Binary  col  1 isions 

4.  Long  correlation  time 

Rp  >  R  >  Rc  Intermediate  region  -  difficult  to  model 

1.  Many  body  problem 

R  >  R  Statistical  broadening  [46]  -  far  wing 

c  1.  Binary  col  1 isions  ~R 

2.  Duration  of  collision  x  ~— 

3.  Short  correlation  time 


j  FW(oj)  describes  the  far  wing  line  shape.  jciu(w)  describes 
the  intermediate  wing  line  shape.  jcwi  p ( t*>)  describes  the  near 
line  center  line  shape.  The  subscriptsLtherefore  represent  the 
frequency  regions  over  which  the  line  shape  function  will  apply. 
This  is  illustrated  in  Figure  12.  The  goal  of  this  section  is 


Figure  12.  The  line  shape  function. 


to  find  the  line  shape  functions  for  near  line  center,  intermediate 
wing,  and  far  wing.  The  near  line  center  and  far  wing  cases  are 
tractable  but  the  intermediate  wing  is  not.  It  will  be  determined 
by  a  mixing  of  the  near  line  center  and  far  wing  line  shape  functions. 

a.  Interruption  broadening 

Table  1  lists  the  conditions  for  interruption  or  impact 
broadening.  These  will  be  applied  to  a  general  correlation 

function,  C  .  (t),  which  describes  the  interactions  between  the 

a  j  D 

absorbing  molecule  and  some  broadener.  Since  the  interaction 
time  is  brief  we  will  reexamine  the  change  in  the  correlation 
function  for  a  small  change  in  time.  Begin  with  Equation  (58) 
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( x+dr) 


usb><Hsb|scaib(T+dT)  |«s1» 


Ro’v’V 


(105) 
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(106) 


SCaib(T+dT)  =  Sc(T-T+dl)Sc(T) 


where 


t  +dx 


S_(x-*x+dx)  =  expl-^x  /  H  .  ( t )  dt 


rn  j  ca-ib' 
o 


and 


t  -  time  at  which  a  collision  occurs, 
o 


Therefore, 


Ca,b(  T+dT)  ’  |Sci,b(  T<T+dT)  |  us^<lst>  |Sca,b<  T*T+dT)  1 ls^\.v.t0 

*«^b|s^Tb^>|^b>«sb|Sca1b(T)|1Lsb»R0.v,to 


(107) 


To  obtain  the  separate  averages  in  the  above  equation  it  is  assumed 
the  events  at  t  and  x+dx  are  independent  [36].  Because  S_,  .(x) 

is  averaged  over  orientation  and  does  not  change  vibrational  and 
electronic  states  for  near  line  center  phenomena,  it  is  independent 
of  the  sum  over  the  quantum  numbers  of  the  broadening  molecule. 

The  general  correlation  function  can  be  written 


Ca,bO>  ■«“* 


SrI.hO)|uSt><Jl! 


Sca,b^> 


«”»„  v  t  -O08) 
O’  ’  0 


Using  this  result  and  Equation  (107)  we  have 


O^^X+dx)  =  Caib(T)^(T-T+dT)>Ro>v>tc 


(109) 


where 


4>(  T*T  +  dx) 


■X 


p  b<us 

b  s 


Sca,b(T"T+dT) 


us 


t><«.sb|sCaib(T-T+dx)  [asS. 


Examine  a  single  collision  event  at  time  t  to  x+dx  or  equivalently 
to  to  tQ+dx  and  the  corresponding  change  in  the  correlation  function 


d  vt> 


-•?  I 


t*?t  c;,b(T) 

I  di-V- 


(no) 


C;  b(x)  indicates  that  no  sum  over  collision  times  is  performed. 

After  evaluating  the  integral  the  following  differential  equation 
is  obtained 


d  Calb(T> 


dt  <J-4>( T+T+dl)^, 


Ro’v  ' 


(111) 


°a,b  =  0-4>(T>x+dT)>R 

0 


-O  bt 

Ca1b<T>  ■  ca,b<0)  e  '  :  Ca,b(0)  =  ’• 


(112) 


Extracting  the  volume  term  in  Equation  (59)  and  writing  the  average 
over  velocity  and  impact  parameter  as 


<  >n  -*■  T7  <  >n 

Ro’v  V  Ro’v 


leads  to 


N.  -n.0,  .  x 

C^b(T)  =  e  b  a1b  . 


(113) 


The  term  nbaa  b  represents  the  magnitude  of  the  collisional  per¬ 
turbations.  For  nboa  b  small  the  absorbing  molecule  will  remain 
correlated  to  its  x=0^state  longer  than  for  nboa  large.  In 

the  frequency  domain  it  is  related  to  the  halfwidth  at  half  maximum 
as  verified  by  the  Fourier  transform  of  the  correlation  function. 
Let  the  halfwidth  be  denoted  by 
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r  1  r  ,  .  i 

aib  aib  aib  b  dlb  b  aib 

r  i 

where  a  and  a  correspond  to  the  real  and  imaginary  parts  of 
nb°a  b'  ^s'n9  Equation  (104)  the  near  line  center  line  shape 
becomes 


,  v 

JcNLC(u))  =  "  7  .  i  )2~~F~T 

+t°a1b) 


where 


(114) 


p  i 

oi  ,  is  the  halfwidth  at  half  maximum  and  a  k  is  the  line  shift. 
a-|D  a-|D 

The  next  step  in  the  analysis  is  to  solve  for  n  k*  First, 

examine 


Sc  (T-*-T+dx)  =  e  1n 
a  i  b 


(115) 


where 


t  +dt 


oo 

n  =  r  I  dt  H  (t)  -  i  f  dt  H  (t)  . 

{  a,b  71  J  ca-,b 

to  1  1 


(116) 


The  extension  of  the  limits  of  integration  is  valid  since  the 
collision  interactions  drop  off  rapidly  for  increasing  intermolecul ar 
distance,  n  can  now  be  written 
alb 

oa  =Zp  b  <l-<usb|eln|usb)<^sb|e'lrlUsl>>R  (117) 

lb  sb  s  1  1  1  o’ 


where  £  p  .=1  has  been  used.  A  perturbation  expansion  on  the 
b  sb 

time  development  operation  can  now  be  carried  out,  that  is 
<£sb|e"ln|£sb>  =  l-i<fcsb|n|£sb>  -  i<£sb|  n2  |£sb>  +  ... 


,  .  1  2  . 

1  “ 1 0j,  -  2  n£  +  *  ’ 


(118) 
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To  second  order 


in  V 


a,  becomes 
alb 


°alb  • 

(119) 

This  is  a  form  similar  to  that  obtained  by  Anderson  [35] .  pta, 
n,  has  been  evaluated  by  Mizushima  [471  using  for  the  collision 
Hamiltonian  the  intermolecul ar  potentials  of  the  general  form 

Hca  (t>  =  ^ — FT7  and  =  (Rn  +  (vt)V/2  and 

calb  R(t)mj  0 


v  R 


mj-1 


r(mj-l)22~mJ 


r(^2 


(120) 


Of  interest  to  this  study  is  the  temperature  dependence 
of  the  halfwidth.  This  can  be  easily  computed  for  the  first  order 
Anderson  theory.  One  notes  that 


1 

n  a  - 


°a]b  ^  ~2  =  v  ^  where 


—  h 

v  is  the  mean  thermal  velocity.  Using  nb  =  the  temperature 

dependence  of  the  halfwidth  is 


(121) 


This  is  the  standard  dependence  used  in  AFGL  tape  [5]  computations. 
It  is  not  however  what  is  sometimes  observed  [48,491.  As  more 
terms  in  the  perturbation  expansion  are  included  the  temperature 
dependence  will  be  modified. 


Two  basic  broadening  mechanisms  can  result  in  the  modulation 
of  the  line  center-frequency,  energy  level  shifts  and  phase  shifts 
of  the  eigenstates.  These  are  analogous  to  frequency  modulation 
and  phase  modulation.  Anderson's  theory  accounts  for  both  effects 
yet  needs  higher  order  terms.  Another  popular  theory  is  the  phase 
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shift  approx  im.it  ion  [47,50)  which  considers  only  phase  modulation 
but  to  all  orders.  To  see  how  the  phase  shift  approximation  comes 
about,  consider  once  again  the  perturbation  expansion  of  Equation 

(HO) 


<lsVin|lsb>.  H(-i<*sb|n|)lsb»  +  ^(ifysVlls^)  +  yy(..)...  . 
Examine 

<«tsb|n°|  est>>  =  V  <dsb |  n | «- '  sb *><«.  ■  sb '  |n|asb> 
rsb' 

=  <esb| n| +  \  <isb|n|l's^'sb|n|is^  . 

<'22> 

For  the  second  term  to  exist  the  operator  n  must  initially  be 
changing  the  eigenstate  to  some  virtual  state  then  back  to  the 
stationary  eigenstate.  This  describes  an  energy  level  shift  or 
frequency  modulation.  The  first  term  in  Equation  (12?)  indicates 
no  change  in  eigenstate  but  may  produce  a  phase  change,  or  phase 
shift.  To  make  the  phase  shift  approximation  simply  ignore  the 
second  term.  Then 


<lsb|e-1n|isb>-  e-\ 


where 


<jtsb|n|nsb>  so  that 


(123) 


o  ,  now  becomes 

tl  1  l) 


O-e 


KhuVn, 


Ro’v 


(124) 


Letting  .(>  -  n  -g  we  may  write 


5? 


Vb  •  )>  b  ,v 

I  ^.b  s  o* 
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(?(m  .1))  2m  . 3 

S  J  '"mb  > 


(  ^=T>C0**5jrr>  "j 


h  s 


(m  1  )fi 


2 

"j'T 


C3 


"j-3 


— 


m, 

r(J)  J 


1 


[<usb|Vj|u*^-<^sb|Vj|tsb>  ?/l"J 


-I 


This  may  ho  written  in  a  more  compact  form  as 


nij-3 

•l  ph «,  T?rar? 


(125) 


J)  S 


Similarly  the  line  shift  term  becomes 


°«ib  =  V’.b^'^X.v 


1  ?sc 


or 


m  .-3 

y  r 

=  L,  »>  J 


s 


*  ta"(VT>  T 


(126) 


For  dipole-dipole  interactions  the  temperature  dependence  of  the 
halfwidth  becomes 


aa,b  tx  T 


(127) 


For  dipole-quadrupole  interactions 


aa,b  • 


(128) 


The  total  near-line-center  line  shape  becomes 
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— 


JcNLC^ 

or 


which  leads  to 


Re  r 


Re  r 


f  i  _a,  -a 

lai-T  3^3 

e  e  e 


alblT 


(129) 


where 


b.  Statistical  broadening 

The  far  wing  and  to  a  lesser  extent  the  line  center  regions 
of  the  true  line  shape  are  affected  by  statistical  broadening  process 
processes.  Strong  collisions  caused  by  close  encounters  dominate 
this  situation.  Because  no  perturbation  expansion  can  be  used, 
the  phase  shift  approximation  will  be  applied.  The  phase  shift 
term,  g,  (see  Fquation  (116))  can  be  written  as 


V, 


dt 


R(t) 


(130) 


The  mgjor  contribution  of  the  phase  term  comes  in  the  region 


t0<T  f51l.  w^re 


2-m . 


V. 

Vn  .  _J  1  / 

n  D  m.-l\ 
vRq  J 


^r(mrl)2  J 

m  2“ 
m. 
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V 

V 


m-  -  even 

J 
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mj  =  3 
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Again  consider  a  general  correlation  function 


or 


CV<T)  V\<"sl’le,n|''^b|''fn|^Ro,v,to] 


N. 

Ca^b<T>  ’ 


%L  0  K<e  4 


1  + 


b>V' 


o’  ’  0 


thus 


(132) 


The  sum  over  tQ  is  converted  ^o  an  integral  which  gives  its  major 
contribution  between  -  ^and  ,  that  is 


<-"v'V-.  <f 


-R. 


or 


O.e^v^K  A  n  'K-Vf. 


(133) 


Evaluating  the  real  and  imaginary  parts  separately  we  have 
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2R, 


<—  C-cos(nu-nj))>|,  - 


(  3/m.  ^  , 

r<-  i:)C0S(?!K>  mj  -  eve" 
J  J 


4it  t  3/mi  3/m.  l 

=  Jh  J(V.  _V  )  J  \ 
m.vfi'  '  ju  js/ 

\J 


which  may  be  written  compactly  as 

•  2R»  .  3/m, 


TT 

1 


<— (l-cos(nu-njl)>R  =  Yjt 


where 


Vju  ■  <usb|Vj|Usb> 
VjS.  -<lsb|Vj|«.sb> 


J  (>£>] 


Similarly  for  the  line  shift  term 
2R 

<—  s,n<Vr'll>>R„,v  = 


mj-3 


(134) 


=  li(')3/l"j(v.  -v.  )3/mj ! 
m,vTr  '  ju  J  \ 


(  3/m.  ,  , 

Cj  mT)sin(|sT)  mj-even 

J  J 


or  in  a  compact  form 
,2R 


3/m, 


TT 


m.=3 

0 


The  general  correlation  function  can  now  be  written 


(135) 
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(136) 


m./3 

Under  the  far  wing  approximation  j  so  that  the  line  shift 

term  can  be  ignored.  The  far  wing  lineJshape  for  the  general 
correlation  function  becomes 


where 


Writing  the  decaying  exponential  function  in  a  power  series,  an 
asymptotic  expansion  can  be  obtained 


jCFW(^  =  n  Re 


y  <V 

hh  S! 


s=0 


_J _ 

l+3s/m . 

( i  Aw)  J 


re*  i7> 

J 


(138) 


The  temperature  dependence  of  the  term  n^.  is 


3/m. 


Vj  (Vj(T))  J  . 


(139) 


The  far  wing  frequency  dependence  for  the  leading  term  in  the 
expansion  is 


jCFW(a,)ot - V37i.  ’ 

(Aw)  J 

The  total  far  wing  line  shape  becomes 


(140) 
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,  i»  T  -A  x3/*a  -Ah  t3/”l  -A.  t3/*2  -A.  t3/™ 

if  o  a  b,  b?  b. 

=  Re  F  e  e  e  e  ...e 


Which,  after  performing  the  Fourier  transform,  becomes 


(141) 
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In  the  very  far  wing  only  the  leading  terms  will  be  necessary. 

It  is  important  to  mention  another  statistical  theory  which 
does  not  use  the  phase  shift  (phase  modulation)  approximation 
by  Fomin  and  Tvorogov  [41].  However  a  perturbation  expansion  only 
to  second  order  is  used.  The  resulting  statistical  line  shape 
contains  the  phase  shift  approximation  formula  multiplied  by  another 
factor.  This  new  term  can  be  interpreted  as  a  level  shift  (fre¬ 
quency  modulation)  correction  factor  to  the  phase  shift  approxi¬ 
mation  result.  Since  energy  level  shifts  are  an  important  part 
of  far  wing  phenomena,  the  correction  term  will  also  be  important. 

A  simplified  expression  for  the  level  shift  correction  term  is 


G.  (m.-l )/m.  2/3 
exp(-  (Av)  J  J)  ) 


(143) 


where  G-apio  •,  p  is  the  reduced  mass  and  w  .  is  a  parameter  which 
decreases  asJj  increases.  Also  G.  does  vary  with  the  vibrational 
quantum  number.  J 


?..  Inhomogeneous  line  shapes 


a.  Doppler  line  shape 


The  random  kinetic  motion  of  molecules  causes  a  Doppler 
shift  of  the  homogeneously  broadened  line  shape.  Molecules  with 
different  velocities  will  experience  different  frequency  shifts 
so  that  a  new  total  line  shape  results,  called  the  inhomogeneous 
Doppler  prof i le. 


The  analysis  of  this  effect  begins  with  Doppler's  equation 


c(v~-v) 


v  =  (c-v)v0  -  v  =  — 


(144) 


and  the  Maxwel 1-Boltzmann  distribution 


WE)  *  e'E/kT  . 


(145) 
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Using  the  kinetic  energy  E  =  mv  and  Equations  (144)  and  (145) 
we  have 


fMB(v-Vo)  =  e 


(v_vo)2 

2kTv20 


(146) 


where  M  is  the  mass  of  the  absorbing  molecule.  This  represents 
the  probability  distribution  of  the  shifted  frequencies  due  to 
velocity  components  of  the  molecules  parallel  to  the  incident 
photon  field.  Enforcing  the  normalizing  condition 


CJO 

\  Jo(v-vo)dv  =  1 
0 


where 


2kTln  2, 


1/2 


is  the  halfwidth  at  half  intensity. 


b.  Voigt  line  shape 

The  Doppler  line  shape  assumes  the  shifted  homogeneous  line 
shapes  are  delta  functions,  which  is  not  strictly  true  but  a  good 
approximation  in  many  cases  such  as  low  total  gas  pressures. 

The  broadened  natural  line  shape  is  always  present  in  any  experiment 
and  in  problems  of  atmospheric  importance,  collision  broadening 
is  present. 

Figure  13  illustrates  the  shifted  collision  broadened  homo- 


Figure  13.  The  Voight  line  shape. 


geneous  line  shape  weighted  by  a  Gaussian  (Doppler)  distribution. 
This  process  can  be  thought  of  as  a  convolution  of  a  homogeneous 
line  shape  through  a  Doppler  profile.  Thus 


Jv(^v0) 


(148) 


60 


L 


Jh(v-Vq)  and  aH  are  the  homogeneous  line  shape  and  its  corresponding 
halfwiath,  respectively.  The  resulting  line  shape  is  called  the 
Voigt  line  shape. 

The  homogeneous  line  shape  that  is  often  used  is  the  Lorentz 
profile.  Substituting  in  Lorentz  and  Doppler  line  shape  the  Voigt 
profile  becomes. 


J'v(v-v0)  = 


°H 

— -j  (In  2) 
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oo 

I 
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where 


v‘  -v 
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E.  The  Absorption  Coefficient  for  H^O  in  N2 


This  section  combines  the  results  of  sections  C  and  D  to 
form  the  absorption  coefficient  particular  to  the  LUO-Np  absorption 
problem.  Initially  the  line  strength  and  line  shape  will  be 
developed  separately. 


The  discussion  of  the  HpO-^  line  shape  must  begin  with 
the  nature  of  the  intermolecuiar  potential  functions.  As  mentioned 
in  Section  B  of  this  chapter  the  intermolecul ar  potentials  are 
not  completely  understood.  This  is  especially  true  concerning 
the  nature  of  strong  collisions  which  form  the  statistical  line 
shape  f 34 ] .  Therefore  a  careful  analysis  concentrating  on  the 
essential  features  of  the  interaction  is  required.  The  general 
form  of  the  intermolecuiar  potential  Hamiltonian  for  H?0-H?0 
collisions,  Hba^,  and  for  N?-H20  collisions,  Hba  ,  becbme, ‘using 

Equation  (14)  through  Equation  (19), 


a  cl 
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Rep2  * 

(151) 


The  subscripts  1  and  2  represent  H~0  and  N2  terms,  respectively. 

Next  the  potential  functions  will  be  considered  for  interruption  and 
statistical  type  collisions. 

The  short  collision  times  characteristic  of  interruption 
broadening  allow  the  approximation  of  Equation  (17)  (i.e., 

<Vp,>  s  Vp.)  to  be  made.  Based  on  the  modeling  success  of  Benedict 
ancl  Kaplan  [52,53]  only  the  leading  terms  in  the  interaction 
potentials  need  be  retained.  Therefore,  for  near  line  center 
phenomena  the  interaction  Hamiltonians  become 


and 


(152) 


(153) 


Collisions  contributing  to  statistical  broadening  may  or 
may  not  require  the  weighted  rotational  average,  <  >r0j--  This 
can  be  determined  by  examining  the  collision  duration  time  and 
time  of  rotation  of  the  molecule.  Based  on  the  average  thermal 
velocities  and  impact  parameters  for  strong  collisions  for  HpO- 
H-0  and  N2-H20  collisions  order  of  magnitude  estimates  for  tne 
col  lision^duration  time  can  be  obtained.  The  result  is 
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10'12 

for  H^O-H^O 

10'13-10"12 

for  N^-H^O 

(154) 


Typical  period  of  rotation  times _^r  ra^ge  between  10"  to  10 

sec  and  for  H~0  range  between  10~  toMO"  sec.  These  numbers 
are  based  on  The  molecular  moments  of  inertia  and  a  representative 
range  of  J-numbers.  Water  vapor  however  can  effectively  rotate 
faster  due  to  its  distortions.  It  is  well  known  that  the  rotational 
motion  of  water  vapor  requires  high  order  perturbation  theory 
[41.  This  necessity  is  a  manifestation  of  the  ease  with  which 
the  water  molecule  will  distort. 


Based  on  the  above  arguments  the  weighted  rotational  average, 
c  >D0t»  is  important  for  hUO-HpO  collisions  and  the  approximation 
of  Equation  (17)  can  be  used  for  N~-Hp0  collisions  in  the  statistical 
limit.  Thus,  for  far  wing  phenomena, The  interaction  Hamiltonians 
become  (retaining  only  leading  terms) 


and 


aa. 


2p^a]  +  2«ie2<r^> 


R3(t) 


Rot 
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Eisenberg  and  Kauzmann  (341  list  the  statistical  average  of 
(^2/R3(t))f  for  rotating  molecules  to  be 


o  4 

^1 _ 

3kTR6( t) 


(157) 


under  the  condition  that 


2Ml 


kTR3 ( t ) 


<< 


1. 
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(158) 


? -  L(x) 

R3(t) 


where  L(x)  the  Langevin  function  and 


x  =  2p‘j/R3(  t) . 


Figure  15 


Figure  15.  Plot  of  the  Langevin  function  L(x). 


is  a  plot  of  L(x).  As  x  increases  above  1  the  curve  becomes  non¬ 
linear  requiring  higher  powers  in  x. 


Because  of  the  complex  nature  of 


2^1 

for  -J—  >  1 
Rot  RJ ( t ) 


and  itSolimited  range  of  importance  (3  to  6  A;  bonding  occurs 
at  2.8  A)  Equation  (157)  will  be  used.  It  is  recognized  that 
temperature  and  frequency  dependent  correction  terms  may  be  required 
in  the  very  far  wing.  Under  this  approximation,  the  H?0-H?0 
interaction  Flamiltonian  becomes  d  d 
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or*  in  a  more  useful  form 
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whore  the  term 


2  2 

e  <rl>i- 


‘1 


2 

represents  a  general  polarizabi 1 ity  for  H?0.  The  quantity  is 

the  or ientat iona 1  polarizabi 1 ity  and  the  remaining  expression  is  the 
electronic  polarizabi 1 ity.  Figure  16  [65]  demonstrates  the  dominance 


Figure  16.  Polarization  of  gases  versus  temperature. 

<t  >r ientat ional  polarizabi 1 ity  over  electronic  polar izabil ity  since 
,.r ientat  ional  polarizabi  1  ity  is  the  only  temperature  dependent 
t.-.m  in  Equation  (160).  Based  on  this,  the  interaction 

*  '  I  f  i  m  t  an  i  s  wr  i  tten 


66 


‘I  3kTR  (t) 


The  near  line  center  and  far  wing  line  shapes  now  follow  from 
Equations  (129)  and  (142), 


J‘«(U)  ■  '  (v-vJW)2 


(162) 


,  .6266  Xa,  .8491  XKA 

i  f  =  I  _ 36  .  _ b4_ 

w  ’ ' "  CT775 


(163) 


where 


■  <  ^>-83[°<^>-1?  Pa  *  pJ 


(164) 


p,  -  absorber  pressure 

a 

p^  -  broadner  pressure 
T  -  temperature 

Xa6  =  FfZh  Pb  (i(V6u"V6lJ)  r(“  co$ty 

sD  s 

(165) 
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Because  of  the  statistical  averaging  over  rotations  the  H-0-H?0 
interaction  potential  matrix  elements  for  strong  collisions  are 
independent  of  angle  and  not  a  function  of  the  J  quantum  number. 
However  Al(v)  will  be  different  for  every  vibrational  band,  identified 
by  the  general  vibrational  quantum  v  representing  both  upper  and 
lower  levels.  For  liL-H-O  collisions  the  statistical  average  was 
not  performed; thus  Az(v,J)  is  a  function  of  the  vibrational  and 
rotational  quantum  numbers  involved.  Since  the  potential  function 
for  near  line  center  and  far  wing  phenomena  are  the  same,  the 
J-dependence  of  A2(v,0)  can  be  expressed  In  terms  of  the  J-dependence 
of  a^.  Using  Equations  (125)  and  (134),  A2(v,J)  becomes 


(167) 


where  .07  represents  a  mean  value  for  a^.  Al(v)  and  A2(v)  will 
be  experimentally  determined  because  of°the  complex  nature  pf 
the  terms  and  uncertainty  of  the  interaction  potentials.  will 
be  obtained  from  the  AFGL  line  listing  (51  and  a  B  of  5  (30?  will 
be  used. 

Equations  (162)  and  (163)  describe  only  parts  of  the  total 
line  shape.  As  section  0  indicates  the  intermediate  wing  region 
is  unknown.  Thus,  the  development  of  a  total  line  shape  must 
take  a  semi-empirical  course.  .The  construction  of  the  total  line 
shape  will  follow  the  guidelines  listed  below. 

1.  The  total  line  shape  is  smooth  everywhere. 

2.  The  total  line  shape  is  normalizable. 

3.  The  near  line  center  and  far  wing  forms  must  dominate 
under  the  appropriate  frequency  limits. 

4.  The  intermediate  wing  region  is  to  some  extent 
a  mixing  of  interruption  and  statistical 
broadening  mechanisms. 
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The  first  step  is  to  make  Equation  (163)  a  bounded  function  for 
all  (v-v  ).  Varanasi  et  al  [56]  have  developed  a  line  shape  which 
accomplishes  this  and  produces  the  statistical  frequency  dependence 
in  the  far  wing.  It  is 


j(v)  .  2  (j?  sin  1)  _  ,  4v=^v0  .  (,68) 

Using  this  result  and  normalizing  the  far  wing  line  shape  to  one 

oo 

(i.e.,  J  JcFW^v^dv=^  Equation  (163)  becomes 


%W 
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(169) 


To  construct  the  total  line  shape  the  concepts  developed  in  the 
introduction  to  Section  0  will  be  followed.  By  choosing  band 
pass  and  band  stop  filters  (centered  at  v  )  with  slow  fall  off, 
mixing  of  the  interruption  and  statistical  line  shapes  respectively 
can  be  achieved  in  the  intermediate  region.  The  total  line  shape 
is  written 


Jc(v)  =  N(JcNLC(v)  P(v)  +  JcFW(  v)  ( 1  - P (  v))  (170) 

where  N  is  the  normalization  constant  determined  by 


J  dv  jc(v)  =  1 
o 

and  p(v)  is  a  band  pass  filter  function  about  v  .  Since  p(v) 
cannot  be  predicted  with  the  present  theory  it  must  be  chosen 
empirically.  To  insure  a  smooth  joining  of  the  interruption  and 
statistical  line  shapes  the  derivative  of  p(v)  at  v=v  and  t  a 
is  required  to  be  zero.  Where  "a"  is  the  frequency  which  represents 
the  cut-off  of  the  window  function  in  the  frequency  domain.  Also 
it  is  desirable  for  p(  v)  to  be  a  well  known  and  easily  integrated 
function.  This  will  be  helpful  in  the  determination  of  N,  the 
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normalization  constant.  Based  on  these  guidelines  p(v)  has  been 
chosen  to  be 


p(v)  =  < 

where  a=5  cm"  . 


+  cos  Av  |  Av |  <5  cm 

0  |Av|>5  cm"^ 

Figure  17  is  a  plot  of  p(v).  As  the  figure 


(171) 


Figure  17.  p(  v)  versus  Au 


illustrates  the  band  pass  filter  function,  p(v),  is  virtually 
unity  a  few  halfwidths  from  line  center.  This  allows  the  inter¬ 
ruption  line  shape  to  dominate  near  line  center  as  it  should. 
Similarly  the  stop  band  filter  function  (l-p(v))  allows  the  statis¬ 
tical  (far  wing)  line  shape  to  dominate  beyond  5  cm"  (50  to  100 
halfwidths  from  line  center).  Thus  a  total  line  shape  has  been 
constructed  which  satisfies  the  four  guidelines  outlined  earlier. 

An  important  correction  to  the  statistical  line  shape  is 
the  level  shift  correction  term  given  by  Equation  (143).  Applying 
this  term  to  Equation  (169)  we  obtain. 
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(172) 


The  A\>  dependence  in  the  exponential  correction  terms  are  slightly 
different  fgom  the  term  derived  by  Fomin  Ul]  (m .=6=> | Av| 
m.=4=>|Av[‘  ).  Both  Av's  are  taken  here  to  the  dne-half  power 
to  simplify  the  normalization  calculations.  The  reduced  mass 
for  Hp0-Hp0  collisions  is  less  than  the  reduced  mass  for  H?0-N? 
collisions.  Also,  from  Equation  (143)  d  c 

“04  >  “06 
and  therefore 

Ga6(v)  <  Gb4<v)  *  073) 

The  "G"  parameters  are  functions  of  the  vibrational  quantum  numbers. 
The  temperature  dependence  of  the  "G"  parameters  is 


(174) 


To  determine  the  line  strength  for  water  vapor,  the  partition 
function  ratios,  S.  ,  E  and  v  in  Equation  (102)  must  be  obtained. 
S^u*  anc*  vo  1)6  *aken  from  AFGL  line  compilation  [5]. 

The  vast  majority  of  the  Hp0  molecules  will  be  in  the  electronic 
ground  state,  therefore  c 


WE1 

*E7 


(175) 


The  vibrational  partition  function  varies  only  slightly  over  the 
temperature  range  of  the  troposphere  and  can  be  assumed  constant. 
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The  rotational  partition  function  for  H~0  is  very  important. 
For  an  asymmetric  top  [57]  the  rotational  partition  function 
ratio  becomes 


c 

QRot 


1.5 


(177) 


The  total  absorption  coefficient  for  H?0-N?  mixtures  can 
be  written  ^  c 


Mv)  =Ik.(v)  =  £  S.u(v)jr(v).  (178) 

i  1  8,u  X,u  c 

The  summation  is  over  all  water  vapor  transitions.  The  line  strength 
in  Equation  (103)  contains  frequency  dependent  terms  which  belong 
to  the  line  shape  in  order  to  satisfy  Equation  (6).  Therefore 
using  Equations  (103),  (162),  (170),  (172),  (175),  (176),  and 
(177)  the  H^O-N^  absorption  coefficient  becomes 


-P(v))  (179) 


where 
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The  norma 1 izat ion  constant,  N,  is  determined  by  the  equation 


CO  oo 

N_1  =J  jCNLC^v^  P(v)dv  +  J  JcfW(W  (1-P(v) )dv  . 


The  solution  of  the  above  equation  as  a  function  of  the  different 
line  shape  parameters  is  presented  in  Appendix  A. 


CHAPTER  III 

EXPERIMENTAL  APPARATUS 


This  chapter  describes  the  experimental  equipment  used  to 
obtain  the  data  of  Chapter  IV.  All  measurements  have  been  taken 
on  one  of  two  long-path  multiple-traversal  absorption  cells  using 
a  CO  laser,  a  C0«  laser  or  a  Fourier  transform  interferometer 
as  probes.  The  combination  of  laser  and  Fourier  transform  spectro¬ 
scopy  techniques  provides  a  more  complete  experimental  environment. 

A.  Frequency  Probes 

1 .  The  CO  laser 

Lasing  action  in  liquid  nitrogen  (LN~)  cooled  CO  lasers 
has  been  observed  from  5pm  to  8pm  corresponding  to  the  1-0  to 
36-35  vibrational  bands  [58].  Due  to  our  interest  in  the  5pm  atmos¬ 
pheric  window,  the  range  of  operation  of  The  Ohio  State  University 
CO  laser  was  optimized  for  operation  on  the  lower  vibrational 
bands. 


The  resonant  cavity  structure  follows  the  design  of  Charles 
Freed  [59].  A  black  granite-invar  structure  supports  the  optics 
and  plasma  tube  and  provides  a  stable  open  loop  optical  system. 

The  four,  1  inch  diameter  Invar  rods,  which  are  thermally  and 
acoustically  shielded,  define  the  length  spacing  between  the  2 
inch  thick  black  granite  end  plates.  Figures  18a  and  18b  show 
a  top  and  side  view  respectively  of  the  Invar  structure.  The 
figures  illustrate  the  different  stages  of  insulation,  acoustic, 
thermal,  and  radiation  shield.  Aluminum  clamps  are  used  to  secure 
the  output  mirror  end  plate  and  Invar  cross  pieces  (used  to  support 
the  plasma  tube  and  liquid  nitrogen  fill  system).  The  grating 
end  plate  is  fastened  to  the  Invar  rods  by  half  inch  stainless 
steel  bolts.  The  output  mirror  endplate  is  9"  x  10-12"  and  the 
grating  end  plate  is  9"  x  12-1/2". 

The  resonant  cavity  optical  system  contains  a  spherical 
output  mirror  of  5  meter  radius  of  curvature  and  a  gold  coated 
grating  mounted  externally  to  the  plasma  tube.  The  spacing  between 
the  mirror  and  grating  is  1.63  meters  producing  an  output  spot 
size  of  2.2  mm  and  a  beam  waist  of  1.87  mm.  The  grating  can  be 
easily  rotated  with  a  micrometer  adjustment  making  the  laser  line- 
selectable.  A  piezoelectric-driven  output  mirror  mount  provides 
the  necessary  adjustments  for  optimum  output  power  and  open  loop 
stabi 1 ity. 
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The  plasma  tube  design  is  based  on  work  by  Djeu  [60].  The 
plasma  is  composed  of  a  N2-He-C0-Xe  mixture  at  a  total  pressure 
of  7  torr.  The  tube  has  nigh  voltage  electrodes  at  each  end  with 
the  ground  electrode  at  the  center  of  the  tube.  Figure  18c  shows 
one  end  of  the  uninsulated  plasma  tube  with  an  electrode.  The 
outer  4"  diameter  pyrex  tube  forms  the  liquid  nitrogen  (LN2) 
reservoir.  The  inner  15  mm  diameter  pyrex  glass  tube  contains 
the  plasma  and  high  voltage  nickel  electrodes.  The  inner  tube 
has  a  glass  bellows  at  each  end  to  ease  the  strain  created  between 
the  inner  and  outer  tubes  when  LN~  is  present.  The  electrodes 
are  1"  long  and  approximately  1"  outer  diameter  and  0.35"  inner 
diameter  with  six  equally  spaced  1/16"  rods  on  the  front  face. 

The  electrodes  not  only  provide  the  high  electric  field  needed 
to  generate  a  plasma,  but  also  serve  to  resist  flow  of  the  CO 
molecules  in  the  active  LN~  cooled  section  to  the  passive  outer 
section.  Thus  ground  state  CO  absorption  is  reduced  to  allow 
operation  of  the  laser  in  the  lower  (2-1  and  1-0)  vibrational 
bands.  The  electrical  connection  to  the  nickel  electrodes  also 
serves  as  a  gas  line  for  the  N2-C0-Xe  mixture  feeding  the  plasma. 
Helium  is  introduced  in  the  warm  passive  section  of  the  tube  and 
serves  as  a  flush  for  unexcited  CO  as  well  as  part  of  the  active 
media.  Figure  18d  shows  the  inlet  ports;  the  top  port  goes  to 
the  high  voltage  electrode  and  the  side  port  is  for  the  helium. 

The  gases  are  then  evacuated  at  the  center  port  near  the  grounded 
electrode.  Also  shown  is  an  end  heater  which  warms  the  exterior 
tubes  to  prevent  moisture  condensation  on  the  CaF?  Brewster  windows. 
The  inlet  lines  from  the  flowmeters  to  the  inlet  ports  on  the 
tube  are  1/16"  inner  diameter;  this  discourages  gas  breakdown 
in  the  inlet  lines  and  makes  the  plasma  more  stable. 

Liquid  nitrogen  CO  lasers  have  the  danger  of  ozone  explosions 
if  the  discharge  is  initiated  with  ozone  present.  Ozone  is  produced 
by  dissociation  of  CO  in  the  plasma,  and  explosions  can  be  a  problem 
when  the  discharge  goes  out  but  the  high  voltage  remains  on. 
Normally,  the  plasma  is  generated  before  LN~  is  introduced  in 
the  reservoir.  For  safe  operation  a  control  circuit  was  constructed 
that  senses  the  state  of  the  plasma  (on  or  off)  and  if  the  plasma 
de-energizes  the  high  voltage  power  supply  and  IN2  fill  control  are 
also  turned  off  automatically.  Figure  19  displays  the  circuit 
schematic  for  this  system.  Since  the  OSU  CO  laser  does  not  use 
02  in  the  active  media  as  other  systems  do,  the  03  formation  problem 
is  greatly  reduced.  Thus,  with  reasonable  care  safe  operation 
of  the  laser  can  be  accomplished. 


One  end  of  the  CO  Laser  plasma  tube  showing  the  electrode 
and  glass  bellow. 


Fig.  18d.  Brewster  window  with  inlet  ports 
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Table  2 

Observed  CO  Laser  Lines 


1 

L  ine 

ID 

v0(cm_1) 

Typical  Power  (mw) 

2-1 

P  ( 9 ) 

2081.258 

2 

P(  10) 

2077.140 

5 

P(ll) 

2072.987 

5 

P(  12) 

2068.802 

7 

3-2 

P(  7) 

2063.225 

1 

P  ( 8 ) 

2059.209 

2 

P(9) 

2055. 159 

10 

P(  10) 

2051.075 

20 

P(H) 

2046.958 

20 

P(12) 

2042.808 

22 

P(  13) 

2038.624 

25 

P(  14) 

2034.4075 

20 

P(  15) 

2030.158 

5 

4-3 

P  ( 7 ) 

2037.124 

5 

P  ( 8 ) 

2033.143 

10 

P(9) 

2029.128 

10 

P(  10) 

2025.080 

10 

P(  1 1) 

2020.998 

40 

P(  12) 

2016.882 

25 

P  ( 1 3 ) 

2012.734 

30 

P(  14) 

2008.552 

25 

P(15) 

2004.337 

25 

P(  16) 

2000.090 

20 

5-4 

P  ( 7 ) 

2011.093 

5 

P  ( 8 ) 

2007.147 

15 

P  ( 9 ) 

2003.167 

15 

P(  10) 

1999.154 

90 

P(  1 1) 

1995.107 

90 

P(  12) 

1991.026 

50 

P(  13) 

1986.913 

40 

P(  14) 

1982.766 

60 

P(  15) 

1978.586 

20 

P(16) 

1974.374 

20 

6-5 

P(6) 

1989.010 

1 

P  ( 7 ) 

1985.133 

6 

P  ( 8 ) 

1981.222 

30 

P  ( 9 ) 

1977.277 

40 

P(  10) 

1973.299 

60 

P(  11) 

1969.287 

60 

P(12) 

1965.242 

40 

P(  13) 

1961.163 

60 

P(14) 

1957.051 

60 

P(  15) 

1952.907 

50 

P(  16) 

1948.729 

50 

P(  17) 

1944.519 

30 

I 
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Table  2  (Continued) 


Line 

ID 

v0(cm  ) 

Typical  Power  (mw) 

7-6 

P(6) 

1963.088 

3 

P(7) 

1959.246 

20 

P(8) 

1955.370 

30 

P(9) 

1951.460 

50 

P(  10) 

1947.517 

60 

P(H) 

1943.540 

50 

P12) 

1939.529 

100 

P(  13) 

1935.486 

60 

P(  14) 

1931.409 

40 

P05) 

1927.299 

15 

8-7 

P(6) 

1937.239 

10 

P(7) 

1933.432 

25 

P(8) 

1929.592 

10 

P(9) 

1925.717  l 

25 

P(  10) 

1921.808 

25 

P(ll) 

1917.866 

25 

P(12) 

1913.891 

25 

P(13) 

1909.883 

20 

P(  14) 

1905.841 

20 

P(  15) 

1901.766 

10 

P(  16) 

197.659 

1 

P(  1 7) 

1893.519 

1 

9-8 

P(6) 

1911.467 

3 

P  ( 7 ) 

1907.695  1 

10 

P  ( 8 ) 

1903.889 

10 

P  ( 9 ) 

1900.049 

25 

P(  10) 

1896.176 

30 

P(  11) 

1892.269 

30 

P(  12) 

1888.328  ! 

25 

P(  13) 

1884.355 

20 

P(  14) 

1880.348 

20 

P(  15) 

1876.309 

20 

P(  16) 

1872.236 

15 

10-9 

P(6) 

1885.771 

4 

P  ( 7 ) 

1882.034 

15 

P  ( 8 ) 

1878.263 

20 

P  ( 9 ) 

1870.620 

40 
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Single  line  operation  has  been  observed  on  8?  CO  transitions 
from  the  ?- 1  vibrational  band  to  the  10-9  vibrational  band.  Research 
efforts  have  emphasized  the  low  vibrational  bands  and  no  attempt 
was  made  to  achieve  output  beyond  the  10-9  band.  It  should  be 
possible  to  extend  the  range  of  the  present  laser  by  two  more 
vibrational  bands.  Outfitting  the  laser  with  different  output 
mirrors  and  gratings  would  extend  the  output  range  to  8  pm.  Table 
2  lists  the  lines  observed  to  date,  along  with  a  nominal  output 
power.  The  CO  laser  presently  operates  with  good  (1%)  short  term 
stability  in  frequency  and  power. 

2.  The  C0o  Laser 

The  CO.,  laser  used  in  the  experiments  has  been  previously 
described  1181.  Briefly,  the  laser  is  a  GTE  Sylvania  model  948 
COp  laser  plasma  tube  with  a  modified  resonant  cavity.  The  laser 
is  also  electronically  stabilized  by  sensing  plasma  impedance 
fluctuations  caused  by  changes  in  the  gain  of  the  system  1611. 

A  feedback  signal  to  the  PZT  cavity  length-controller  maintains 
optimum  mirror  separation.  Lasing  action  on  23  transitions  in 
the  00  l-lO^O  (10pm)  vibrational  band  and  20  transitions  in  the 
00°l-02°0  (9pm)  vibrational  band  has  been  observed..  This  effectively 
covers  a  frequency  region  from  930  cm'1  to  1085  cm*  ,  ideally 
suited  for  10pm  window  studies. 

3.  The  Fourier  Transform  Spectrometer 

The  monochromatic  outputs  of  the  laser  sources  prevent  the 
observation  of  the  total  spectral  response  of  a  molecule.  The 
addition  of  a  Nicolet  7199  Fourier  transform  spectrometer  expands 
the  experimental  capability  of  this  study.  The  continuous  spectral 
coverage  from  800  cm  to  5000  cm  with  a  resolution  of  0.05 
cm  (0.035  cm  for  spectra  of  limited  bandwidth)  provided  by 
the  instrument  allows  the  close  examination  of  local  absorption 
line  and  near  band  wing  features  of  absorbing  molecules.  The 
frequency  response  of  the  LN-  cooled  HgCdTe  detector  determines 
the  lower  limit  of  the  spectral  range. 

The  heart  of  the  spectrometer  is  a  Michelson  interferometer 
with  a  fixed  mirror  in  one  arm  and  a  moving  mirror  in  the  other. 

The  moving  mirror  intensity  modulates  a  broadband  collimated  beam 
and  produces  an  interferogram,  I ( ) ,  which  can  be  thought  of  as  a 
sum  of  single  monochromatic  interference  patterns. 


B( \>)  cos(?tt\kS) 


(184) 


1(6) 


/ 


B(v)  cos(2nv6)  dy 


where 

B(u) 

amplitude  component  at  frequency  y  in  cm' 

,*=2L 

optical  retardation  length 

L 

distance  traveled  by  moving  mirror. 

To  obtain  the  spectrum  B(  j)  a  Fourier  transform  of  !(*)  is  taken 


B(  y)=2 


/ 


I(<5)  cos(2irv6)d6 


o 


(185) 


Modern  computers  and  the  fast  Fourier  transform  (FFT)  algorithm 
by  Cooley  and  Tukey  f621  have  made  Fourier  transform  spectroscopy 
a  practical  tool.  However  the  digitization  of  the  interf erogram 
and  finite  optical  retardation  lengths  limit  the  resulting  spectrum 
in  resolution  and  bandwidth.  To  perform  the  FFT  and  other  data 
manipulations  the  system  includes  a  Nicolet  1180  computer.  The 
computer  handles  up  to  one-half  million  point  (complex)  transforms. 

B.  White  Absorption  Cells 

Two  White-type  absorption  cells  have  been  used  in  the  experi¬ 
ments,  a  15.24  meter  base  path  cell  constructed  in  the  early  sixties 
and  a  10.785  meter  base  path  cell  just  recently  completed.  Ex¬ 
perience  with  the  older  absorption  cell  has  prompted  many  improvements 
in  the  10.785  meter  cell.  A  brief  description  of  these  cells 
fol lows. 

1 .  The  15.24  meter  cell 

This  cell  was  first  described  by  Long  f63l  and  more  recently 
by  Peterson  [181.  The  cell  is  a  steel  pipe,  0.607  m  in  diameter, 
and  16.15  m  long.  The  walls  of  the  cell  were  honed  smooth  during 
construction  in  order  to  reduce  water  vapor  adsorption  effects. 

The  mirrors  are  aluminum-coated  pyrex  with  a  measured  reflectivity 
of  97.1%  at  5  gm  and  radii  of  curvature  of  15.24  meters  which 
determines  the  mirror  separation.  The  entrance  and  exit  windows 
are  currently  BaF~  half-degree  wedges.  A  six  inch  diffusion  pump 
and  100  cfm  mechanical  pump  are  used  for  cell  evacuation. 


fhe  long  path  lengths  (1-1.5  km)  required  for  accurate  measure¬ 
ments  of  the  weak  water  vapor  absorption  at  C0?  and  certain  low- 
vibrational  CO  laser  frequencies  prompted  a  careful  study  of  the 
15. ?4  m  cell  stability.  Several  problems  were  identified  and 
are  discussed  below. 

The  optical  alignment  of  the  cell  was  sensitive  to  changes 
in  the  internal  pressure.  The  output  beam  would  shift  upwards 
by  3  mm  at  a  1  km  path  length  as  the  cell's  pressure  was  changed 
from  vacuum  to  one  atmosphere.  Nonuniformities  of  the  windows, 
the  mirrors,  and  the  detector  then  resulted  in  fluctuations  in 
the  output  beam's  signal  level  resulting  in  a  high  noise  level 
on  the  data.  To  solve  this  problem  the  background  measurement 
(usually  taken  with  the  cell  evacuated)  was  taken  with  the  cell 
at  the  same  pressure  as  the  sample  measurements.  This  procedure 
has  the  additional  benefit  of  ratioing  out  absorption  effects 
of  contaminants  in  the  buffer  gas  if  the  same  gas  supply  is  used 
throughout  the  experiment. 

Vibrations  of  the  cell  caused  by  vacuum  pump  and  building 
vibrations  (i.e.,  compressors  switching  on  and  off  etc.)  also 
changed  the  mirror  alignment  in  the  cell.  Because  typical  measure¬ 
ments  take  an  entire  day  it  was  difficult  to  complete  successfully 
a  long-path  experiment.  By  v ibrat ional ly  isolating  all  pumps 
and  the  mixing  fans  inside  the  cell  the  vibration  problem  was 
greatly  reduced.  To  best  avoid  building  noise  the  experiments 
were  performed  during  off  hours  whenever  possible. 

Another  problem  which  had  a  major  effect  on  the  experiment 
was  nonuniform  reflectivity  of  the  mirror  surfaces.  Oil  deposits 
and  particulates  collected  on  the  mirrors  over  several  years  in 
a  random  fashion.  Thus  any  change  in  the  cell  alignment  had  a 
dramatic  effect  on  the  output  signal  level.  The  problem  was  solved 
by  recoating  the  mirrors  and  by  using  a  100  mw  Argon  laser  to 
maintain  alignment  throughout  the  experiment. 

Figure  20  illustrates  the  experimental  arrangement  for  the 
CO  and  C0o  laser  measurements.  Irises  1  and  2  were  used  to  verify 
the  alignment  of  the  ir  lasers  against  the  Argon  alignment  laser. 

A  BaFp  beamsplitter  was  used  to  direct  part  of  the  laser  energy 
to  a  reference  detector.  Mirrors  mounted  on  kinematic  mounts 
are  labeled  KM.  The  focal  lengths  of  spherical  mirrors  are  listed 
next  to  the  mirrors. 

2 .  The  10,785  m  cell 

The  10.8  meter  cell  which  will  now  be  described  was  designed 
by  Professor  Edward  K.  Damon  [641.  Or.  Robert  J.  Nordstrom  conceived 
the  unique  optical  system  which  couples  the  lasers  and  interferometer 
to  the  optical  system  of  the  absorption  cell. 
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The  10.735  m  cell  features  many  design  improvements  over 
the  15.34  m  cell.  The  internal  cell  mirrors  are  mounted  externally 
through  vacuum  bellows  in  order  to  reduce  the  pressure  dependent 
alignment  problem  of  the  older  cell  and  to  maintain  a  constant 
mirror  separation  over  a  wide  temperature  range  of  the  cell. 

The  mirrors  are  microprocessor  controlled  by  stepper  motors  allowing 
great  ease  in  changing  or  correcting  the  cell  alignment.  They 
are  gold  coated  with  a  reflectivity  of  93.6%  at  10.6|#ii,  and  have 
a  radius  of  curvature  of  10.785  m.  The  improved  optical  design 
of  the  10.785  m  cell  features  longer  path  lengths,  improved  spatial 
stability,  and  reduced  aberration  effects  [65,66)  over  the  standard 
White  type  design  [671.  Figure  21  illustrates  the  spot  arrangement 
on  the  field  mirror  of  the  new  cell.  Spot  number  5  would  typically 
be  the  output  spot  on  a  White  type  system.  Instead  it  is  fed 
back  into  the  optical  system  generating  two  additional  rows  of 
spots.  Path  lengths  of  1.8  km  with  laser  sources  and  1  km  with 
a  broadband  source  have  been  achieved  on  the  new  cell.  The  cell 
and  optical  tables  are  rigidly  coupled  and  vibrationaly  isolated 
from  the  building.  The  new  cell  is  constructed  of  304  stainless 
steel  with  a  smooth  interior  finish  to  reduce  adsorption  effects 
and  to  allow  a  wider  variety  of  gases  to  be  studied  (e.g.  ozone). 
Buffer  gases  are  introduced  uniformly  along  the  cell.  This  method 
has  greatly  reduced  gas  mixing  times  in  the  new  cell  over  the 
mixing  time  in  the  older  15.24  m  cell.  The  two  cells  have  similar 
vacuum  systems  except  the  10.735  m  cell  has  a  Key  Vacuum  Product 
coaxial  foreline  trap  model  number  CFT-300  to  reduce  oil  back- 
streaming  of  the  mechancal  Dump.  The  new  cell  is  also  temperature 
controlled  from  -60°C  to  60nC  with  a  length  uniformity  of  rC. 

This  in  fact  was  the  principal  reason  for  designing  the  new  cell. 

Figure  22  illustrates  the  total  experimental  system  of  the 
10.785  m  cell  facility.  The  enclosures  labeled  source,  FTIR, 
transfer  optics,  and  detectors  represent  sealed  boxes  which  house 
the  indicated  equipment.  The  labels  S  and  W  in  the  optical  path 
of  the  blackbody  source  represent  the  source  image  and  the  FTIR 
beamspl itter  image  to  be  imaged  on  the  field  and  focusing  mirrors 
of  the  absorption  cell,  respectively.  The  sliding  mirror  in  the 
detector  box  can  be  inserted  into  the  cell  output  optical  path 
to  test  the  alignment  of  the  optical  system.  The  notation  used 
in  Figure  2.2  is  the  same  as  Figure  20. 

C .  Experimental  Cons  i derat  ions 

White  type  absorption  cells  represent  an  important  part 
of  the  experimental  apparatus  presently  used  for  the  measurement 
of  weak  absorption  coefficients.  Also,  the  spectrophone,  another 
sensitive  experimental  tool,  requires  calibration  which  is  usually 
accomplished  by  White  cell  data.  Because  of  their  importance 
a  detailed  discussion  of  the  nature  of  White  cell  experiments 
fol lows. 


White  cells  measure  the  transm i ttance  through  a  medium  as 
opposed  to  a  spectrophone  which  directly  measures  the  absorption 
effect.  This  can  give  the  spectrophone  an  experimental  advantage 
for  transmi ttances  above  90%  [68].  The  transmittance  is  the  ratio 
of  the  radiation  intensity  attenuated  by  a  medium  to  the  unattenu¬ 
ated  radiation  intensity.  Thus, 


T 


where 


rB(n)7TRB 

T 

1  s  ( n ) 


IB(n) 


086) 


transmittance 

attenuated  intensity  level  due  to  the  sample 

reference  detector  signal  level  during  sample 
measurement 

unattenuated  intensity  level  i background) 

reference  detector  signal  level  during 
background  measurement. 

number  of  traversals. 


IRc  andlg„  scale  I  (n)  and  I B ( n )  respectively  to  account  for  the 
different  power  levels  of  the  source.  The  sample  intensity,  I  (n), 
can  be  further  expressed  as 


rs(n)  =  I0SRr(x,y 

-kn£ 

’PH20^Ks(t)e  (187) 

where  R(x,y,PH^Q) 

-  the  mirror  reflectivity  as  a  function 
of  the  mirror  surface  and  water  vapor 
pressure 

!os 

-  intensity  level  at  input  to  cell  during 
sample  measurement 

Ks(t) 

-  time  dependent  intensity  fluctuations 
caused  by  windows,  amplifiers,  detectors, 
mirror  drift  etc. 

r 

-  number  of  reflections  inside  cell 

2. 

0 

-  base  path  length  of  cell. 

To  eliminate  the  time  dependent  fluctuations  a  time  average  over 
a  number  of  measurements  is  taken.  Two  experimental  techniques 
can  be  used  to  accomplish  this  averaging  process. 
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Since  the  mirror  reflectivity  is  a  function  of  the  position 
on  the  mirror  surface  and  cell  vibrations  can  cause  mirror  drift 
which  misaligns  the  cell  and  changes  the  position  of  the  spots 
on  the  field  mirror,  intensity  variations  in  the  output  beam  can 
result.  Path  differencing  [69 ]  averages  over  the  mirror  surface 
by  changing  the  path  length  while  keeping  the  sample  fixed  and 
averages  over  window,  detector,  etc.  noise  by  taking  many  measure¬ 
ments  at  each  path  length.  Another  approach  used  at  The  Ohio 
State  University  maintains  the  cell  alignment  throughout  an  experi¬ 
ment  by  using  optical  lasers  to  account  for  the  drift  of  the  cell 
mirrors.  By  greatly  reducing  the  mirror  drift  problem  no  average 
over  the  mirror  surface  is  required.  Figures  20  and  22  show  the 
alignment  lasers  and  irises  1  and  2  which  are  used  to  verify  align¬ 
ment  of  the  cell  before  each  measurement  sequence.  Also  the  output 
spot  from  the  cell  is  monitored  to  test  the  cell  mirror  alignment 
before  each  measurement  sequence.  The  average  over  the  other 
noise  sources  is  accomplished  by  averaging  400  data  points  with 
a  45  millisecond  interval  taken  on  a  computer-controlled  data 
acquisition  system  for  a  fixed  sample. 


The  time 

The  background 

^B^X 

Using  Equation 


averaged  intensity  is  expressed  as 

_  .  .  -knj, 

!os  R  ^x,y,pH20^  <Ks(t)\  e 

intensity  signal  is  similarly  written,  ( k  =0 ) 

!oB  R  (x,y,pH20=0^  ‘ 

(186)  the  time  averaged  transmittance  becomes 


(188) 


(189) 


R>(X,y,PH20}  <Ks(t))t  -knn0 

Rr(x,y,PH^0=0) 


where 


1. 


(190) 
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I’ropor  time  averaging  requires 


<Ml)>t  ■  <KBmX  . 


Therefore, 


Rr(  x»y»P(.|^()) 

Rr(x»y*Va0) 


-knC 

(  ) 
e 
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To  obtain  the  true  absorption  coefficient,  k,  (see  Equation  (4)) 
the  mirror  reflectivity  must  be  the  same  in  vacuum  (or  dry  gaseous 
mixture)  as  it  is  with  water  vapor  present  at  any  pressure.  This 
point  has  proven  to  be  a  major  obstacle  in  obtaining  the  true 
absorption  coefficient.  Water  vapor  is  adsorbed  by  the  mirrors 
forming  a  very  thin  film  (>rl  micron  [70').  The  thickness  of  the 
film  will  depend  on  mirror  temperature  and  pressure  of  water  vapor 
in  the  cell.  Therefore  the  effect  on  mirror  reflectivity  also 
depends  on  temperature  and  pressure.  The  reduction  of  mirror 
reflectivity  by  water  vapor  will  also  be  frequency  dependent. 
Interference  effects  become  more  important  as  the  wavelength  approaches 
the  film  thickness  (i.e.,  R  decreases  as  v  inc  eases).  Absorption 
and  refractive  index  effects  exhibit  an  irregular  frequency  dependence 
since  they  are  related  to  resonance  phenomena  (71).  Burch  et 
al  '14'  reported  in  1967  that  water  vapor  adsorption  did  not  cause 
major  difficulties  at  room  temperature  with  10  torr  water  vapor 
and  36  reflections.  However,  Burcii  has  recently  reported  [721 
that  the  4  pm  region  has  major  water  vapor  adsorption  problems. 

Further  the  10.735  m  cell  at  OSU  has  110  to  150  reflections  and 
typically  has  more  than  10  torr  water  vapor  for  long  path  experi¬ 
ments.  Thus,  water  vapor  adsorption  effects  are  important  in 
particular  at  4  ^m  and  to  a  lesser  extent  at  10  pm. 
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■Tivr  A  IV 

ANALYSIS  OF  RESULTS 


A  considerable  amount  of  infrared  window  water  vapor  absorption 
data  has  been  collected  over  the  past  ten  years  at  The  Ohio  State 
University  and  elsewhere.  Measurements  of  room  temperature  nitrogen- 
broadened  water  vapor  absorption  made  during  this  study  and  collected 
during  experimental  programs  at  other  laboratories  will  be  presented. 
The  total  line  shape  expression  developed  in  Chapter  II  will  then 
be  used  in  conjunction  with  the  AFGL  line  listing  [5]  to  model 
the  measured  data.  The  results  of  the  modeling  demonstrate  the 
importance  of  far  wing  line  shape  contributions  to  the  water  vapor 
continuum  absorption. 

A.  Data  Base 


Table  3  is  a  compendium  of  nitrogen-broadened  water  vapor 
window  absorption  paramerers  at  room  temperature  (296  K).  The 
general  form  for  the  absorption  coefficient  which  is  used  to  present 
the  data  [ 1 4 J  is 


Mv)  =  (CN  PN)  +  Cspa)  (km' '  1  ,  (192) 

where  C  is  the  self-broadening  coefficient  for  water  va,or, 

is  the  nitrogen  broadening  coefficient,  pM  is  the  pressure  ot 

in, 

nitrogen  in  torr,  p  is  the  water  vapor  partial  pressure  in  torr, 

R  is  the  ideal  gas  constant,  and  T  is  the  temperature  in  Kelvin. 

The  equation  can  be  conveniently  rewritten  to  obtain 


k(  v) 


CN 

ITT  Pi 


(PN?  +  Bpa} 


(193) 


where  B  =  is  the  dimensionless  self-broadening  coefficient. 

lN  C 

^N 

Table  3  lists  ^y  and  B  at  many  different  frequencies.  One  exception 

to  the  formula  (Equation  (193))  must  be  made  for  1957.050  cm~^ 
since  the  CO  laser  line  at  this  frequency  is  within  the  halfwidth 
of  a  local  water  absorption  line.  In  this  case  the  formula  becomes 


k(v)  =  py  pa(p 


N, 


Bpa). 
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Table  3 


Nitrogen  Broadened,  Water  Vapor 
Window  Absorption  Parameters 


v( cm-  1 ) 

CN/RTxl05(— 

torr^ 

337.9 

656.7 

366.0 

165.2 

389.0 

124.9 

411.0 

60.1 

433.7 

68.7 

439.0 

68.7 

448.8 

50.2 

475.1 

67.4 

482.6 

71.2 

498.0 

23.7 

531.6 

14.4 

559.2 

8.67 

579.0 

23.5 

597.0 

11.5 

611.4 

5.24 

629.0 

5.80 

934.894 

0.36 

936.804 

0.32 

938.688 

0.498 

940.548 

0.571 

942.383 

0.759 

944. 194 

0.257 

945.980 

0.164 

951.192 

0.308 

952.881 

0.331 

969.140 

0.041 

970.547 

1.761 

971.930 

1.091 

973.288 

0.837 

974.622 

0.903 

977.214 

0.768 

980.913 

0.163 

1048.661 

0.555 

1050.441 

0.730 

1052.196 

0.674 

1053.924 

0.588 

lu73.278 

0.401 

1077.303 

0.583 

1079.852 

0.167 

1081.087 

0.564 

1082.296 

0.546 

1854.933 

42.9 

B 

Laser  Line 

ID 

_ 

10.5 

a 

25.4 

a 

26.9 

a 

45.0 

a 

35.6 

a 

35.0 

a 

43.6 

a 

27.4 

a 

23.6 

a 

55.2 

a 

70.8 

a 

86.6 

a 

30.7 

a 

56.2 

a 

95.9 

a 

84.6 

a 

330.0 

P ( 30 )  ,b 

343.6 

P ( 28 )  ,b 

195.7 

P(26)  ,b 

171.5 

P(24)t,b 

114.8 

P( 22)  ,b 

340.1 

P(20)+,b 

549.5 

P(  18)  ,b 

264.6 

P(12)  ,b 

251.2 

P( 10)  ,b 

2133.4 

R( 10)*, b 

35.1 

R(12)  ,b 

107.3 

R(14)*,b 

95.2 

R( 16)*, b 

80.6 

R( 18)  ,b 

98.0 

R(22)  ,b 

467.3 

R (28)  ,b 

120.9 

P(  18)  l",b 

86.9 

PC  16 )  f\b 

99.0 

P(  14) t,b 

138.6 

P(12)+,b 

213.2 

R(12)*,b 

113.2 

R ( 18)  ,b 

414.9 

R(22)*,b 

131.2 

R(24)  ,b 

105.6 

R(26)+,b 

7.6 

10-9  P( 14)  , c 
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Table  3  (rontinuod) 


v ( cm" ' ) 

,  .  -1 

CN /RTxUV’t  -•  iy) 
torr‘ 

8 

Laser 

12 

.  i  ne 

1 880.348 

26.7 

16.  1 

9-8  P ( 14) 

,c 

NOS.  840 

38.9 

9.9 

8-7  P(  14 ) 

,c 

19? 7. 299 

23.4 

10.4 

7-6  P( 15) 

.c 

1931.409 

17.2 

8.2 

7-6  P( 14) 

,c 

1948.729 

19.7 

13.3 

6-5  P( 16) 

.c 

1952.907 

10.5 

15.2 

6-5  P (  16) 

,c 

1957.050 

44.7 

3.8 

6-5  P ( 14) 

,c 

1970. 1 29 

16.8 

13.6 

5-4  P ( 1 7 ) 

.c 

1974.374 

4.88 

26.6 

5-4  P ( 16) 

,c 

1978.586 

4.46 

21.5 

5-4  P( 15) 

.c 

2003. 165 

5.04 

7.25 

5-4  P(9) 

,c 

2004.337 

5.21 

4.4 

4-3  P ( 15) 

,c 

2012.733 

6.65 

5.5 

4-3  f* ( 13) 

,c 

2026.079 

4.95 

5.2 

4-3  P ( 10) 

,c 

2029.128 

2.20 

8.75 

4-3  P(9) 

,c 

2038.625 

5.34 

6.5 

3-2  P( 13) 

.C 

2471.243 

0.26 

54.3 

3-2  P(ll) 

.8 

2496 . 720 

0.25 

44.6 

3-2  P( 10) 

.d 

26:’ 1.769 

0.233 

37.3 

3-2  P( 9) 

,d 

2546.373 

0.221 

32. 1 

3-2  P(8) 

,d 

2663.951 

0.216 

31.4 

2-1  P( 1 1 ) 

.d 

2570.522 

0.210 

29.7 

3-2  P( 7) 

,d 

2580.095 

0.204 

29.7 

2-1  P( 10) 

,d 

2594.197 

0.  198 

30.3 

3-2  P ( 6 ) 

,d 

2605.806 

0.191 

32. 1 

2-1  P ( 9 ) 

,d 

2617.386 

0.186  1 

34.5 

3-2  P(5) 

.d 

2631.066 

0.179 

38.6 

2-1  P(8) 

,d 

2640.075 

0.175 

42.0 

3-2  P ( 4 ) 

,d 

2655.863 

0.168 

50.2 

2-1  P( 7) 

,d 

2665.218 

0.  165 

55.2 

1-0  P ( 10) 

,d 

2680.178 

0. 157 

67.  1 

2-1  P(6) 

,d 

2691.608 

0.152 

78.1 

1-0  P(9) 

,d 

2703.998 

0.146 

90.9 

2-1  P(5) 

,d 

2717.538 

0.141 

107.5 

1-0  P(8) 

,d 

2727.308 

0.137 

121.9 

2-1  P ( 4 ) 

.  d 

2742.997 

0.132 

144.9 

1-0  P( 7 ) 

,  d 

2750.093 

0.128 

158.7 

2-1  P(3) 

,d 

2767.968 

0.123 

192.3 

1-0  P ( 6 ) 

.d 

2792.434 

0.116 

250.0 

1-0  P(5) 

.d 

2816.380 

0. 1 10 

312.5 

1-0  P ( 4 ) 

,d 

2839.791 

0.102 

400.0 

1-0  P ( 3 ) 

,d 

2862.646 

0.099 

476.2 

1-0  P(2) 

,  d 

3434.9994 

13.3 

5.00 

2-1  P ( 8 ) 

,e 
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I  .ill  1 1 •  !  (Cont  inued) 


v(cm~ 1 ) 

C  /R  1x10s (-k-~=) 
torr‘ 

B 

Laser 

L  ine 

ID 

3483.652? 

69.0 

8.20 

2-1  P( 7) 

,e 

3531.  1/47 

476.  7 

3.91 

2-1  P  6) 

,e 

3577.500? 

532.8 

6.30 

7-1  P ( 5 ) 

,e 

3644.1454 

429.4 

5 . 00 

1-0  P( 7) 

,e 

3693.4226 

1542.3 

4 . 34 

1-0  P ( 6 ) 

,e 

3788.7753 

203.4 

8.9? 

1-0  P ( 4 ) 

,e 

Sources  of  experimental  data. 

a.  Gryvnak  and  Burch  17  1) 

h.  Peterson,  Thomas,  Nordstrom,  Damon  and  Long  [701 

r .  lorn),  Damon,  Nordstrom,  Peterson,  Thomas  and  Sherman  (731 

d.  Watkins,  White,  Bower  and  Sojka  [341 

e.  Watkins,  Spell  icy.  White,  Sojka,  Bower  |7hJ 

The  values  from  337.9  to  679.0  can  *  are  taken  from  Burch 
[ 7 3 1  wfio  used  a  (jrat.  ing  spectrometer  and  a  White-t^ee  absorption 
cell.  The  values  from  934.894  cm  to  1082.796  cm  "  represent 
t fit'  00°l-10°0  and  00°l-0?l\)  CO,  laser  line  measurement s  taken 
at  The  Ohio  State  University  1701.  The  asterisk  in  the  laser 
lint'  identification  column  indicates  the  measurement  using  The 
Ohio  State  University  differential  spectrophone  [201;  the  remainder 
of  The  Ohio  State  Uni  ersity  measurements  were  determined  using  the 
lb. 24  m  White  cell  l70l.  The  complete  set  of  spectrophone  measure¬ 
ments  show  excellent  agreement  with  the  White  cell  data  [181. 

Further  verification  of  the  listed  values  designated  by  a  dagger 
(f)  in  the  laser  line  identification  column  has  been  achieved 
through  a  recently  completed  remeasurement.  program  using  the  newly 
completed  10.786  m  cell  which  was  described, in  the  last  chapter. 

The  data  from  1854.933  cm  to  7038.675  cm’  are  CO  laser  measure¬ 
ments  also  taken  at  The  Ohio  State  University  (23,74)  on  the 
15.74  in  cell.  The  parameters  of  the  last  seven  CO  laser  lines 
were  obtained  with  the  Cl^  laser  described^ in  Chapter  111.  The 
values  from  7471.743  cm  to  7867.646  cm"  are  DF  laser  measurements 
performed  by  the  Atmospheric  Sciences  Laboratory  (74,751 ._jThe 
II F  laser  measurements  from  3434.9994  cm’  to  3788.7753  cm  were 
also  performed  by  the  Atmospheric  Sciences  Laboratory  [761. 
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The  magnitude  of  absorption  at  a  particular  frequency  can 
be  inferred  from  the  magnitude  of  CL/RT  and  B.  The  water  vapor 
partial  pressure  dependence  of  the  absorption  coefficient  can 
be  determined  from  the  magnitude  of  B.  For  B  small  (^5-10)  the 
k  versus  p  dependence  is  nearly  linear  and  for  B  large  the  dependence 
is  nearly  parabolic.  The  table  clearly  shows  the  characteri sties 
of  the  different  frequency  regions.  The  interpretation  and  modeling 
of  this  data  is  the  topic  of  the  next  section. 

To  complement  the  monochromatic  frequency  listing  of  Table 
3,  broadband  spectral  data  of  pure  water  vapor  have  also  been 
taken  using  the  Nicolet  FTIR  system  (see  Chapter  II j).  The  portions 
of  this  spectrum  in  the  window  regions  from  850  cm  to  1282  cm'  , 

1800  cm*  to  2054  cm"  ,  and  2500  cm”  to  2862  cm'  are  presented 
in  Appendix  B.  This  FTIR  spectrum  is  for  a  sample  of  12  torr 
pure  water  vapor  at  296  K  and  323.5  meters  path  length.  The  local 
line  structure  in  the  windows  is  clearly  indicated.  The  resolution 
is  0.05  cm"  . 

B.  Application  of  the  Total  Line  Shape 

This  section  contains  three  parts.  The  first  examines  the 
far  wing  "A"  and  "G"  parameters  and  the  normalization  coefficient 
of  the  total  line  shape  developed  in  Chapter  II.  The  determination 
of  the  "A"  and  "G"  parameters  allows  the  application  of  the  total 
line  shape  formulas  to  the  water  vapor  continuum  problem.  This 
application  becomes  the  topic  of  the  next  two  parts. 

Part  two  presents  a  comparison  of  experimental  and  theoretical 
absorption  coefficients  versus  frequency  and  FLO  partial  pressure. 

The  agreement  is  good,  from  which  it  is  concluded  that  far  wings 
can  explain  the  observed  FLO  continuum  absorption.  The  window 
regions  studied  are  the  4  (jm,  5  vm,  10  pm,  17  pm  and  millimeter. 

The  general  success  of  the  model  in  predicting  absorption  in  all 
the  FLO  infrared  windows  is  further  evidence  of  the  correctness 
of  the  total  line  shape  approach.  Shortcomings  in  the  model  are 
also  recognized  and  discussed. 

Part  three  presents  the  observed  and  modeled  temperature 
dependence  of  the  absorption  coefficient.  The  far  wing  expression 
for  the  total  line  shape  features  a  negative  temperature  dependence 
which  is  required  if  agreement  with  experiment  is  to  be  obtained. 
However  the  rate  of  the  negative  temperature  dependence  is  less 
than  observed  experimentally.  This  is  explained  by  recognizing 
the  inadequacies  of  the  temperature  dependence  of  the  intermolecular 
potentials. 


t 


95 


1 .  Par.neters  of  the  total  line  shape 


The  "A"  and  "G"  far  wing  parameters  of  the  total  line  shape 
(Equation  (182))  will  be  determined  from  experimental  data.  To 
avoid  any  complicating  effects  (i.e.,  dimer  absorption)  the  experi¬ 
mental  data  are  taken  from  "near-band"  regions  where  far  wing 
contributions  from  the  band  dominate.  Furthermore,  this 
allows  the  determination  of  the  far  wing  contribution  in  the  12-8  urn 
and  4  t im  window  regions.  Figure  1  shows  15  pm,  5  pm  and  3  pm 
to  be  near-band  regions  for  H^O  absorption.  The  15  .m  and  longer 
wavelength  data  by  Burch  [72]  are  used  to  obtain  the  line  shape 
parameters  for  the  rotational  band.  The  5  pm  CO  laser  results 
are  used  to  determine  the  parameters  for  the  \>2  fundamental  band. 

The  3  pm  HF  laser  results  are  used  to  obtain  tne  parameters  for 
the  pj  and  fundamental  bands  (see  Table  3). 

A  computer  program  applying  the  total  line  shape  equations 
and  employing  the  AFGL  absorption  line  listing  for  water  vapor 
was  used  to  model  the  experimental  data.  The  parameters  were 
determined  by  trial  and  error  until  a  good  representation  of  the 
near  band  experimental  data  was  obtained.  The  success  of  this 
technique  removed  the  need  for  more  sophisticated  fitting  routines. 
Further,  considering  the  limitations  of  the  model,  a  more  accurate 
determination  of  the  "A"  and  "G"  parameters  would  not  be  meaningful. 
The  "A"  and  "G"  parameters  for  the  main  H,0  bands  are  listed  in 
Table  4.  c 


Table  4 


The  Far  Wing  Parameters 


A1 (Rotational )  =  2.35 

G, (Rotational )  =  0.065 

A2(Rotational )  =0.18 

Gb(Rotational )  =  0.115 

A1(v2)  =  1 

.73 

Ga(v2)  =  0.07 

A2(v2)  = 

.153 

Gb(v2)  =  0.084 

AUx^.v-j) 

=  1.95 

Ga( V1 ’ v3 )  =  *045 

A2(vj,v3) 

=  0.165 

Gb( Vi ,v3)  =  .088 

Using  the  equations  of  Appendix  A  and  the  far  wing  parameters 
of  Table  4,  the  normalization  coefficient,  N,  can  be  evaluated 
for  the  different  water  vapor  bands.  For  water  vapor-nitrogen 
mixtures  under  general  tropospheric  conditions  of  temperature 
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and  pressure  and  for  y^lOO  cin"  the  normal  izat  ion  coefficient 
is  equal  to  one  within  one  percent.  At  less  than  100  cin  the 
normal izat ion  coefficient  is  less  than  one  (i.e.,  v=  18.6  cm" r 
N=0.95).  Since  most  of  the  water  vapor  spectral  lines  are  above 
100  cm  ,  N  is  generally  equal  to  one.  The  total  line  shape  contains 
a  normalized  near- 1 ine-center  line  shape  (Equation  (181)).  When 
multiplied  by  the  p-function  it  represents  98%  to  96%  of  the  area 
under  the  total  line  shape,  very  close  to  what  the  near- 1 ine-center 
area  would  be  if  no  p-function  were  used.  Therefore,  the  total 
line  shape  as  constructed  for  the  different  bands  does  not  modify 
the  emphasis  of  the  near  line  center  line  shapes  commonly  used 
(the  Lorentz  line  shape  for  v^lOOO  cm  and  the  modified  Van  Vieck- 
Weisskopf  (451  line  shape  for  \«1000  cm"1).  This  is  an  important 
point  because  previous  modeling  attempts  in  the  5  t<n  region  [74] 
have  used  super-Lorentz  line  shapes  which  enhance  the  wings  of 
the  local  lines  while  ignoring  the  far  wings  of  strong  lines  in 
the  heart  of  the  fundamental  band.  To  maintain  normal  izat  ion 
the  absorption  at  Tine  center  must  decrease  by  a  significant  amount 
(10-15%).  This  distorts  the  true  total  line  shape  since  the  Lorentz 
line  shape  does  model  near  line  center  phenomena  (Chapter  II 
Section  D).  The  line  shape  developed  in  this  study  avoids  this 
problem. 

To  reduce  the  long  computer  time  of  an  all  line  calculation 
requires  two  conditions  to  be  met.  One  is  a  bound  (frequency  range 
in  which  absorption  lines  are  considered  about  a  central  frequency) 
of  1500  cm  .  Larger  bounds  produce  no  significant  change.  Note 
that  previous  worker^  with  the  Lorentz  line  shape  have  used  a 
bound  of  only  20  cm  which  apparently  was  selected  to  minimize 
computer  time  rather  than  from  any  physical  principles.  Also 
our  program  performs  a  test  on  the  absorption  line  intensity  to 
eliminate  the  many  weak  lines  which  do  not  have  appreciable  far 
wing  absorption.  The  test  criteria  is 
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where  S?  and  are  read  directly  from  the  AFGL  listing  and  A\> 
is  the  frequency  difference  between  the  absorption  line  and  the 
frequency  location  for  which  the  absorption  coefficient  is  to 
be  computed.  If  Equation  (195)  is  satisfied  then  the  absorption 
line  is  ignored. 


n 
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frequency  and  pressure  dependence 


The  frequency  dependence  of  the  statistical  line  shape  far 
wing  (Equation  (1 82))  is  illustrated  in  Figure  23  for  the  rotational, 

Vp  and  (v,,v,)bands  of  nitrogen  broadened  H«0.  Comparison  to 
the  Lorentz  line  shape  is  also  made.  The  figure  shows  the  statis¬ 
tical  fa^  wing  to  be  greater  than  the  Lorentz  far  wing  until  around 
1000  cm  (this  intersection  point  decreases  with  pressure)  from 
line  center  where  the  exponential  factor  in  the  statistical  line 
shape  causes  a  drop  below  the  Lorentz  far  wing.  Thus  the  total 
line  shape  (Equation  (179))  can  be  technically  termed  sub-Lorentzian. 

The  results  of  the  application  of  the  total  line  shape  determined 
for  each  band  are  presented  in  Figures  24  thru  30.  R *  c ? 1 T  that 
figures  24,  25,  28  and  30  represent  the  near  band  data  from  which 
the  "A”  and  ”G"  parameters  were  determined.  Figures  26,  27  and 
29  represent  the  window  regions  and  show  remarkable  agreement 
with  experimental  data.  This  is  a  strong  indication  that  far 
wings  of  the  water  vapor  fundamental  bands  play  an  important  role 
in  continuum  absorption.  Discrepancies  in  the  data  comparison 
can  be  attributed  to  experimental  error,  inaccurate  parameters 
on  the  AFGL  [5]  listing,  and  breakdown  of  the  line  shape  model. 

Perhaps  the  most  serious  breakdown  of  the  model  occurs  in  the 
10  in  window  region,  which  is  dominated  by  the  rotational  band. 
Collisj^nal  interactions  resulting  in  line  distortions  600  to 
800  cm  from  line  center  are  not  small  perturbations  to  the  pure 
rotational  motion  of  the  molecule  which  will  typically  have  energy 
levels  from  0  to  a  few  hundred  wavenumbers  (see  Chapter  2,  Section  B). 
Thus  the  model  of  the  far  wings  for  the  rotational  band  is  not 
strictly  valid.  However,  the  results  in  the  10  iim  region  are 
surprisingly  good  and  perhaps  this  is  not  a  serious  problem. 

Also,  the  inability  to  achieve  a  reasonably  good  fit  to  the  HF 
data  makes  the  "A"  and  "G"  parameters  for  the  bands  question¬ 

able. 

Random  errors  in  the  absorption  line  shape  parameters  (line 
position,  line  strength,  halfwidth,  etc.)  tend  to  average  out 
in  the  far  wing  continuum  region  because  hundreds  of  absorption 
lines  are  contributing  simultaneously.  However  local  line  contri¬ 
butions  do  not  have  this  property.  Appendix  B  contains  spectra 
from  850  to  2862  cm"  of  the  window  regions  studied.  The  local 
line  structure  is  clearly  indicated  and  a  correlation  can  be  found 
between  strong  local  line  contributions  and  data  discrepancies 
in  many  cases.  The  problem  of  properly  characterizing  nearby 
absorption  lines  relative  to  some  frequency  of  concern  may  not 
only  pertain  to  the  standard  line  shape  parameters  (i.e.,  line 
strength,  halfwidth  and  line  position).  Since  line  positions 
are  usually  obtained  from  low  pressure  samples  of  pure  water  vapor 
modeling  attempts  at  a  total  pressure  of  one  atmosphere  may  require 
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Figure  24.  Comparison  of  exoerimental  and  theoretical  frequency 
dependence  of  the  absorption  coefficient 
from  300  cm"  to  440  cm"  . 
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Figure  26.  Comparison  of  experimental  and  theoretical 
dependence  of  the  absorption  coefficient 
from  930  cm"  to  990  cm. 
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Figure  27.  Comparison  of  experimental  and  theoretical  frequency 
dependence  of  the  absorption  coefficient 
from  1040  cm'  to  1090  cm. 
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of  the  absorption  coefficient  from  1850  cm  to  2050  cm’ 
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the  consideration  of  pressure  shifts.  Tin's  effect  is  described 
in  Chapter  II,  Section  D  (see  Equation  (126)).  The  pressure  or 
line  shift  is  similar  to  the  halfwidth  (pressure  and  temperature 
dependence)  except  smaller  in  magnitude.  In  fact  line  shifts 
are  usually  ignored  because  it  is  a  small  effect.  However  even 
a  small  change  in  line  position  of  a  local  absorption  line  can 
have  significant  effects  on  the  absorption  characteristics  at 
a  monochromatic  frequency.  This  is  especially  true  if  the  probe 
frequency  is  within  two  halfwidths  of  the  line  center  position. 

This  condition  is  quite  often  satisfied  in  the  near  band  regions. 

Figure  29  illustrates  the  4  jin  window  region.  The  calculated 
data  points  are  of  two  types.  One  set  represents  the  total  amount 
of  absorption,  local  line  and  far  wing  contributions.  The  other 
set  represents  only  far  wing  absorption.  Comparison  to  experimental 
measurements  by  Burch  (151  is  made  showing  reasonable  agreement 
with  the  far  wing  calculations.  Since  Burch  made  measurements 
where  local  line  contributions  were  minimal,  only  far  wing  contri¬ 
butions  are  represented.  Table  5  compares  experimental  OF  laser 
measurements  by  Mills  177)  to  calculations  using  the  total  line 
shape  of  this  study.  The  measurements  include  local  line  contri¬ 
butions  and  the  modeling  agrees  quite  well  in  every  case.  The 
OF  laser  measurements  (24)  listed  in  Table  3  are  generally  higher 
than  measurements  made  by  Burch  and  Mills.  A  comparison  between 
the  measurements  of  Burch  and  ASL  is  made  in  Figure  3.  As  emphasized 
in  Chapter  III,  White  cell  measurements  in  this  region  are  very 
difficult  and  discrepancies  between  experimenters  are  expected. 
However,  the  measurements  by  Burch  and  Mills  and  the  far  wing 
calculations  show  good  consistency. 


Table  5 

Comparison  of  OF  laser  measurements  by  Mills 
to  the  total  line  shape  calculations. 


T  =  296  K  p  =  14.3  Torr  Pw  =  745.7  Torr 

o  N 


V 

Calculated  (km-^) 

Mills  (km' ^ ) 

2546.373 

.027 

.031 

2570.522 

.032 

.035 

2594.197 

.037 

.057 

2631.066 

.033 

.030 

2655. 063 

.103 

.103 

2600. 179 

.067 

.075 
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The  water  vapor  partial  pressure  dependence  of  the  absorption 
coefficient  is  characterized  by  the  magnitude  of  B  in  Equation 
(19.1).  For  B  small  (5-10)  the  linear  term  will  dominate  and  for 
B  large  (100-200)  the  quadratic  term  will  dominate.  The  curves 
of  the  absorption  coefficient  versus  water  vapor  partial  pressure 
can  vary  from  linear  to  parabolic.  Using  the  far  wing  formula 
of  the  total  line  shape  (Equation  (182))  an  expression  for  B  can 
be  obtained 


_  .3198  A1 ( v) 

'  7TOT  ‘A?(Y;j) 


.5 


exp(  (Gb4-Ga6)»/Av)  (Ax») 


1/4 


(196) 


Figure  31  shows  a  plot  of  B  versus  A\>  for  the  different  bands. 


Figure  31.  Frequency  dependence  of  the  self -broadening 
coefficient  for  the  different  bands. 


Notice  that  the  choice  of  the  "A"  and  "G"  parameters  has  produced 
larger  B  values  for  the  rotational  and  bands  than  the 

band.  This  feature  is  consistent  with  wnat^is  observed  in  liquid 
water  spectra. 
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Since  liquids  are  composed  of  closely  spaced  molecules 
(relative  to  a  gas)  the  nature  of  the  collisions  will  always  be 
in  the  statistical  domain.  Thus  liquid  spectra  provide  an  oppor¬ 
tunity  to  examine  general  features  of  the  statistical  line  shape 
as  well  as  the  physics  of  close  collisions.  Figure  32  is  a  plot 
of  a  liquid  water  spectrum  from  800  to  4800  cm"  .  (t  shows  the 
corresponding  fundamental  bands  of  the  gas  phase,  the  v?  bending 
mode  and  v,,\u  stretching  modes.  Two  other  vibrat  ionaTbands 
also  appear.  JThe  librational  band,  around  800  cm"  on  the 
plot  is  caused  by  the  hindered  rotational  motion  of  the  H~0  molecules 
[341.  The  association  band,  v. ,  at  2050  cm"  is  caused  by  the 
hindered  translational  -  rotational  motion  of  the  H~0  molecules  [341. 
The  perturbation  of  the  rotational  motion  of  H?0  is  represented 
by  the  average  over  orientations  <  >R  .  as  discussed  in  Chapter 
11,  Section  B.  It  is  clear  from  the  H~0  liquid  spectrum  that 
this  is  an  important  consideration.  Further  the  band  strength 
of  the  band  is  significantly  larger  than  the  band  as 

opposed  to  the  gas  phase  where  they  are  more  nearly  equal.  This 
suggests  that  strong  collisions  perturb  the  stretching  motion 
of  H~0  more  than  the  bending  motion.  This  fact  explains  the  different 
magnitudes  of  B  for  the  v?  and  vi,v3  bands  in  the  statistical 
wing.  The  rotational  band  is  also  greatly  perturbed  by  the  strong 
collisions  and  therefore  has  large  B  values  similar  to  the  v-i,v, 
bands.  1  J 


Table  3  lists  experimental  B  values  at  many  frequencies. 

Notice  that  as  the  frequency  increases  away  from  the  rotational 
and  bands  B  increases  as  one  would  expect.  At  337.9  cm'1  B 
is  10.5  then  increases  to  200-300  in  the  10  vm  window.  The  fluctuation 
in  B  indicated  by  Table  3  to  small  values  at  unexpected  frequencies 
can  be  explained  by  local  line  effects.  Since  near  line  center 
B  is  5  this  will  decrease  the  effective  B  that  is  observed.  Notice 
also  that  in  the  4  ^  window  region  B  is  lower  on  the  side  closest 
to  the  \u  band.  This  is  consistent  with  the  observation  that 
B  will  be  smaller  for  band  than  the  v-),^  bands. 

As  an  example  the  effective  value  of  B  determined  from  least 
squares  fitting  calculated  points  to  Equation  (193)  for  the  00°  1- 
10  0  P(20)  laser  line  position  is  123.  This  number  represents 
a  weighted  average  of  B  values  from  many  absorption  lines.  The 
major  contributing  lines  from  the  rotational  band  are  centered 
around  200  cm  and  therefore  from  Figure  31  the  B  value  for  these 
lines  is  approximately  between  170  and  200.  The  resulting  B  value 
is  less  than  this  because  of  contributions  from  the  vh  band  (see 
Figure  31)  and  local  lines  which  have  smaller  B  values. 
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The  success  of  the  model  in  predicting  the  proper  water 
vapor  partial  pressure  dependence  of  the  absorption  coefficient 
is  demonstrated  by  Figures  33  thru  46.  The  plots  are  of  nitrogen- 
broadened  water  vapor  absorption  at  ?96l K  and  760  Torr  total  pressure. 
The  solid  curve  is  the  theoretical  prediction  and  the  boxes  are 
experimental  data  points.  The  frequency  of  the  absorption  is 
labeled  on  each  plot.  The  frequency  locations  selected  for  the 
plots  have  good  magnitude  agreement  between  theory  and  experiment 
emphas ir inc  the  models  ability  to  predict  the  pressure  dependence. 

The  millimeter  and  submillimeter  spectral  regions  in  the 
troposphere  also  exhibit  continuum  absorption  [78].  The  total 
line  shape  model  should  also  work  in  this  region.  Figures  47 
thru  49  represent  window  absorption  at  millimeter  frequencies. 

The  solid  curves  are  theoretical  calculations  and  the  squares 
are  experimental  data  points  taken  from  Bast  in  [78].  The  lack 
of  experimental  data  in  this  region  makes  it  impossible  to  char¬ 
acterize  this  region  more  completely.  The  model  predicts  the 
absorption  magnitude  fairly  well  at  7.667  cin  but  shows  increasing 
relative  absorption  as  the  frequency  decreases.  This  general 
trend  can  be  caused  by  improper  frequency  dependence  of  the  rotational 
band  red  wing  as  it  goes  to  0  cm”  .  A  strong  candidate  for  this 
problem  is  the  symmetric  structure  of  the  statistical  line  shape 
(Equation  (169)).  The  asymmetry  of  distortion  of  a  low  lying 
rotational  energy  level  in  a  vibrational  potential  well  should 
produce  an  asymmetric  line  shape.  The  total  line  shape  model 
used  in  this  study  does  not  account  for  this. 

3 .  Temperature  dependence 

The  temperature  dependence  of  continuum  absorption  has  been 
a  difficult  problem.  As  pointed  out  in  Figure  14  close  collisions 
which  generate  the  far  wings  violate  the  approximation  made  on 
the  interaction  potential.  Thus  it  can  be  expected  that  the  tempera¬ 
ture  dependence  in  the  far  wing  will  not  be  totally  valid.  The 
problem  is  further  complicated  for  the  rotational  band  far  wings 
because  of  the  breakdown  in  the  perturbation  expansion  as  discussed 
previously.  A  correction  to  the  self-broadening  coefficient  of 
the  rotational  band  lines  is  used  in  the  model.  However  no  modifi¬ 
cation  to  the  temperature  dependent  factors  in  the  total  line 
shape  equations  governing  the  vibrational  bands  has  been  used. 

Table  5  lists  the  C../RT  and  B  experimental  coefficients 
[.73]  for  the  water  vapor  absorption  coefficient  at  430  K  from 
430.0  cm’  to  6B9.0  cm'  .  The  values  of  B  in  this  frequency  region 
are  on  the  average  a  factor  of  6  smaller  than  the  corresponding 
R  values  in  Table  3.  The  temperature  dependence  of  B  indicated 
in  Equation  (196)  does  not  account  for  this.  Since  this  region 
is  dominated  by  far  wing  behavior  and  it  is  recognized  that  the 
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Figure  35.  Comparison  of  experimental  and  theoretical 
H^O  partial  pressure  dependence  at  611.4  cm"  . 
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Table  5 


Nitrogen  Broadened  Water  Vapor  Window 
Absorption  Parameter  at  430  K. 


V 

Cn/RT  X  105  (km'1  J 
Torr^ 

B 

430.0 

103.1 

6.9 

433.7 

190.5 

7.0 

440.1 

111.6 

7.8 

448.8 

127.6 

7.9 

465.4 

87.4 

7.8 

475.1 

103.1 

7.6 

482.6 

104.3 

6.6 

498.0 

32.5 

10.9 

531.6 

18.3 

12.3 

559.2 

13.0 

13.1 

579.0 

40.7 

7.0 

597.0 

19.2 

10.3 

629.0 

8.5 

13.0 

far  wing  temperature  dependence  is  not  accurately  represented, 
a  correction  to  Equation  (196)  was  performed.  By  changing  the 

power  of  (— r-)  from  0.5  to  4.5  in  Equation  (196)  a  substantial 


improvement  in  modeling  performance  was  achieved.  That  particular 
term  comes  from  the  temperature  dependence  of  the  interaction 
potential  (see  Equation  (157)).  Figure  50  illustrates  the  agreement 
between  the  calculated  and  experimental  data.  Figures  51  and 
52  show  the  HpO  partial  pressure  dependence  at  430  K.  Excellent 
agreement  with  experimental  measurements  is  demonstrated. 


The  application  of  the  total  line  shape  to  the  temperature 
dependence  at  the  00°1-10°0  P(20)  laser  line  is  illustrated  in 
Figure  53.  The  ordinate  is  C  as  defined  in  Equation  (192). 

The  two  labeled  solid  lines  represent  calculations  using  the  total 
line  shape  and  the  Lorentz  line  shape  with  all  the  absorption 
lines  included  from  the  rotational  and  Vp  bands  (i.e.,  no  bound). 

The  Lorentz  line  shape  curve  shows  a  slight  positive  temperature 
dependence  and  is  considerably  lower  than  the  experimental  data. 

The  total  line  shape  calculation  does  not  show  the  rapid  temperature 
dependence  of  the  experimental  data  but  the  computed  C  does  decrease 
with  increasing  temperature.  The  experiment  1  data  points  are 
taken  from  Peterson  (OSU  data)  fl8]  and  Arefev  and  Dianov-Klokov 
(Soviet  data)  [211.  High  temperature  measurements  by  Montgomery  [221 
at  1203  cm"  exhibit  an  interesting  effect.  A  negative  temperature 


129 


53*. 6  cm  1 


Figure  52.  Comparison  of  experimental  and  theoretical  H.,0 
partial  pressure  dependence  at  531.6  cm"T  at  430  K. 
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Figure  53.  Comparison  of  experimental  and  theoretical 
temperature  dependence  of  C$  at  944.1945  cm"^. 


dependence  of  C  is  observed  from  290K  to  400K,  then  a  positive 
dependence  appears  above  400K.  A  minimum  occurs  around  400K. 

Figure  54  shows  a  comparison  between  Montgomery's  data  and  calculations 
using  the  total  line  shape  formula  developed  in  Chapter  II.  The 
general  trends  of  the  temperature  dependence  are  modeled.  However 
the  magnitude  of  the  temperature  dependence  is  not  closely  modeled. 

Burch. [151  indicates  a  decrease  in  C  from  296K  to  430K 
at  2600  cm  by  a  factor  of  3.  The  total  Tine  shape  model  predicts 
a  factor  of  1.2  decrease. 

Continuum  modeling  by  far  wings  of  Lorentz  pr  Van  Vleck- 
Weisshopf  line  shapes  indicate  a  positive  temperature  dependence 
in  C  .  This  fact  has  been  used  to  discredit  far  wing  models  in 
general.  The  important  result  of  this  section  is  that  a  properly 
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CHAPTER  V 
CON Cl  US  IONS 


Attenuation  of  infrared  radiation  in  the  troposphere  is 
dominated  by  water  vapor  absorption.  Most  past  work  on  the  spectro¬ 
scopy  of  water  vapor  has  been  on  the  analysis  of  the  rotational 
and  vibrational  bands.  As  a  result  a  thorough  listing  of  line 
positions,  line  strengths  and  halfwidths  exists  today.  Because 
of  the  weak  "continuum"  absorption  in  the  water  vapor  windows 
centered  at  10  and  4  microns,  efforts  to  measure  and  model  the 
continuum  have  not  been  as  successful  as  efforts  in  the  analysis 
of  Hie  bands.  Yet  a  precise  understanding  of  continuum  absorption 
is  important  for  long  path  energy  transmission.  For  this  reason 
the  study  of  the  water  vapor  windows  lias  had  a  long  and  speculative 
h i story . 

1  lie  purpose  of  this  study  is  to  demonstrate  the  importance 
of  far  wing  phenomena  in  characterizing  H90  continuum  absorption. 

A  total  line  shape  for  water  v apor-ni trogbn  interactions  valid 
under  tropospheric  conditions  is  derived.  The  model  is  tested 
using  a  set  of  experimental  room  temperature  H.,0  continuum  measure¬ 
ments  of  high  quality.  Using  this  data  base  parameters  of  the 
far  wing  component  of  the  total  line  shape  are  determined  from 
near  band  experimental  data.  Grating  spectrometer  measurements 
from  300  to  6B0  cm  are  used  to  determine  unknown  far  wing  param¬ 
eters  of  the  pure  rotationa)  band  of  H.,0.  CO  and  HF  laser  measure¬ 
ments  taken  in  the  5  micron  and  3  mierhn  regions  are  used  to 
determine  the  far  wing  parameters  of  the  \>0  and  fundamental 

bands,  respectively. 

The  total  line  shape  model  is  applied  to  the  millimeter, 

10  pm,  and  4  pm  absorption  windows  with  encouraging  success. 

A  significant  increase  in  the  self-broadening  ability  of  H?0  over 
N,  is  predicted  in  the  far  wing.  This  allows  the  proper  modeling 
of  the  absorption  coefficient  versus  H.,0  partial  pressure  dependence 
in  all  window  regions.  A  negative  temperature  dependence  is  predic¬ 
ted  by  the  model  in  the  continuum.  The  observed  rate  of  the  tempera¬ 
ture  decrease  is  not  predicted  by  the  model,  however  the  failings 
are  related  to  approx imat ions  made  on  the  interaction  potentials 
and  perturbation  expansion  of  the  Hamiltonians.  Although  the 
total  line  shape  model  has  limitations,  it  does  demonstrate  the 
importance  of  considering  far  wings  of  absorption  lines  in  continuum 
absorpt ion. 


Included  in  the  study  are  spectra  of  hLO  in  the  10  pm, 

|  in  and  4  |in  window  regions  obtained  with  a( Fourier  transform 
spectrometer.  The  resolution  is  0.05  cm  .  The  spectral  plots 
•ire  presented  in  Appendix  B. 

f  it  i  •  ••fforts  to  improve  the  total  1  ini'  shape  model  should 
concentrate  on  the  statistical  line  shape  part.  The  phase  shift 
approximation  which  is  used  does  not  adequately  describe  the  nature 
of  the  statistical  far  wings.  The  level  shift  correction  term 
is  an  important  factor  in  the  far  wing  equations.  A  more  general 
solution  including  phase  and  level  shifts  of  the  stationary  states 
should  be  attempted.  Also  future  investigations  of  the  molecular 
interactions  important  in  far  wing  physics  must  solve  the  problem 
of  strong  interaction  energies  (greater  than  kT). 

The  theoretical  development  in  Chapter  II  presents  a  complete 
characterization  of  the  molecular  absorption  problem.  However 
approximations  made  during  the  solution  making  the  problem  tractable 
have  caused  shortcomings  in  the  model  as  previously  described. 

Ihe  important  contribution  of  Chapter  II  is  the  unified  development 
of  line  strength  and  line  shape  in  a  general  setting.  This  has 
allowed  the  description  of  the  absorption  coefficient  from  millimeter 
to  the  near  infrared. 
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APPENOIX  A 

NORMALIZATION  INTEGRALS 


The  solution  of  the  normal ization  constant  (Equation  (183)) 
is  presented  in  this  Appendix.  This  requires  the  integration 
of  the  total  line  shape  (Equation  (170))  and  making  the  area  under¬ 
neath  the  line  shape  equal  to  one.  The  realization  of  this  require¬ 
ment  resulted  in  the  evaluation  of  many  difficult  integrals. 

Whenever  possible  an  exact  solution  was  obtained.  However  in 
many  cases  numerical  and  series  expansion  approaches  were  used. 

To  simplify  the  computation  of  the  normalization  constant 
for  every  absorption  line  the  total  line  shape  formula  has  been 
written  for  two  cases  \>  <1000  cm  and  \>o>1000  cm  (see  Equations 
(181)  and  (182)).  Also°the  simpler  expression  for  the  case 
v  >1000  cm  reduced  the  computer  time  of  the  program  used  for 
calculating  the  absorption  coefficient. 

A  summary  of  the  resulting  expressions  is  presented  in  the 
fol lowing  equations. 
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APPENDIX  B 

Ho0  SPECTRA  IN  THE  WINDOW  REGIONS 


Appendix  B  presents  H„0  spectra  in  the,window  reqions  from 
HSO  cm  _ j  to  128?  cm'1,  1800W  to  2054  cm'1  and  2500' cm'7  to 
2862  cm'  .  The  plots  represent  a  12  torr  water  vapor  sample  at 
280  K  taken  on  a  Nicolet.7199  Fourier  transform  spectrometer  with 
a  resolution  of  0.05  cm  .  The  ordinate  of  the  following  plots 
is  thepihsorpt ion  coefficient  in  km“  and  the  abscissa  is  frequency 
in  cm'  .  The  measurement  was  made  on  the  10.785  m  cell  with  a 
path  length  of  323.6  m. 
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Figure  59.  H,0  spectra  from  398  cm  to  912  cm 
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Figure  76.  H20  spectra  from  1100  cm  1  to  1114  cm 


LABORATORY  SPECTRUM 


CD 


□s*z.  os  *  s  os*e  os*  i  os*-'- 

33NtJgyosad 


171 


Figure  77.  spectra  from  1112  cm'1  to  1126  cm 
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Figure  78.  i^O  spectra  from  1124  cm  to  1138  cm 
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Figure  79.  rLO  spectra  from  1136  cm  to  1150  cm 
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Figure  31.  Ho0  spectra  from  1160  cm  to  1174  cm 
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Figure  83.  H?0  spectra  from  1134  cm  to  1198  cm 
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Figure  35.  ’i90  spectra  from  1203  cn  to  1222  cm 
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Figure  89.  H2O  spectra  from  1256  cm-i  to  1270  cm 
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Figure  91.  ILO  spectra  from  1309  cm  to  1314  cm 
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Figure  101.  H2O  spectra  from  1920  cm'1  to  1934  cm 
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Figure  105.  H2O  spectra  from  1963  cm  to  1982  cm 
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Figure  107.  i^O  spectra  from  1992  cm  to  2006  cm 
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Figure  111.  HO  spectra  from  2040  cm-1  to  2054  cm 
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Figure  112.  H£ 0  spectra  from  2500  cm-1  to  2514  cm'1. 
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Figure  116.  H,0  spectra  from  2548  cm"1  to  2562  cm 


560  2562  2564  2566  2568  2570 

WAVENUMBERS 

Figure  117.  H20  spectra  from  2560  an'1  to  2574  cm" 
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Figure  119.  H?0  spectra  from  2596  cm  to  2610  cm' 
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Figure  121.  H-0  spectra  from  2620  cm  to  2634  cm 
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Figure  122.  H,0  spectra  from  2632  cm  to  2646  cm 
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LRBORRTORY  SPECTRUM 


spectra  from  2716  cm*  to  2730 
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Figure  133.  HO  spectra  from  2764  cm  to  2778  cm 


Figure  134.  H?0 
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Figure  135.  HO  spectra  from  2793  cm  to  2802  cm 
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Figure  136.  H?0  spectra  from  2300  cm-i  to  2314 
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Figure  140.  H,0  spectra  from  2848  cm"1  to  2862  cm" 


